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PIPING  INELASTIC  FRACTURE  MECHANICS  ANALYSIS 

1.0  SUMMARY 

This  report  summarizes  the  results  and  conclusions  of  Task  1  and  2  of  the  study  on  'Piping  Ine¬ 
lastic  Fracture  Mechanics  Analysis",  Contract  Number  (NRC-03-79-116).  In  these  tasks,  available 
experimental  data  and  the  analytical  methods  for  predicting  rupture  of  LWR  piping  have  been  assem¬ 
bled  and  assessed.  The  analytical  techniques  investigated  can  sbe  catalogued  into  three  major  groups. 
These  are  gross  structural  response  analysis,  semiempirical  methods,'  and  the  J-controlled  growth 
approach. 

The  gross  structural  response  is  computed  using  numerical  techniques  when  dynamic  behavior  of 
both  the  piping  and  the  fluid  is  considered.  The  LWR  piping  stability  is  modeled  by  allowing  a  preas¬ 
signed  crack  to  propagate  under  a  certain  toughness  criterion  such  as  Klc,  K,d  or  maximum  strain  at  the 
crack  tips.  In  this  analysis,  the  accuracy  of  the  results  depends  very  heavily  on  the  crack  growth  cri¬ 
terion  about  which  little  information  is  known.  When  brittle  failure  is  observed,  K,c  or  Ktd  may  be 
adequate.  For  a  pipe  that  fails  plastically  or  in  a  mixed  brittle/ductile  failure  mode,  well  established 
fracture  criteria  are  not  yet  available.  In  addition,  a  single  computer  run  frequently  provides  numerical 
results  which  may  not  be  extrapolated  to  other  cases  where  computing  costs  prevent  additional  runs  to 
be  made. 

Depending  upon  the  piping  materials  and  their  service  temperature,  a  semiempirical  approach  has 
been  proposed  to  study  elastic/plastic  fracture  of  LWR  piping.  A  pseudo-toughness  parameter,  Kc,  has 
been  derived  from  classical  fracture  mechanics  theory  and  modified  to  describe  the  toughness  of  the 
piping  made  of  low  to  medium  toughness  materials.  Unfortunately,  poor  correlation  was  found 

between  the  pseudo-toughness,  Ac,  and  experimental  results.  For  high-toughness  materials,  a  flow 
Manuscript  submitted  April  30,  1980. 


1 


CHANG,  NAKAGAK1,  GRIFFIS  AND  MASUMURA 


-4 


stress  theory  has  been  proposed.  In  this  case,  the  theory  characterizes  the  piping  stability  by  its  'imit 
load  carrying  capability  and  contains  no  fracture  parameter  of  any  type.  By  properly  adjusting  the  flow 
stress  value,  experimental  results  from  tests  on  piping  made  of  high  toughness  materials  suggest  that 
the  flow  stress  theory  may  be  adequate.  However,  upon  further  evaluation,  it  is  revealed  that  structural 
stability  is  influenced  by  both  crack  tip  conditions  as  well  as  structural  geometry.  Therefore,  more 
conclusive  analytical  or  experimental  evidence  is  needed  before  any  conclusion  can  be  drawn  regarding 
the  adequacy  and  the  applicability  of  the  semiempirical  methods  in  LWR  rupture  prediction. 

One  of  the  most  promising  and  most  rigorous  predictive  methods  to  date  is  the  J-control  growth. 
J/c  has  been  accepted  as  an  elastic/plastic  fracture  initiation  criterion  and  the  existence  of  a  HRR 
(Hutchinson-Rice-Rosengren)  held  has  been  identified  as  the  necessary  condition  of  a  J -dominated 
stress  field.  Once  a  crack  starts  to  propagate,  however,  large  deformations  and  unloading  near  the  crack 
tip  region  are  expected;  these  factors  are  not  accounted  for  in  the  original  J-theory  formation  where 
infinitesimal  deformation  and  deformation  theory  of  plasticity  (no  unloading)  are  assumed.  Neverthe¬ 
less,  extensive  research  effort  has  been  directed  to  extend  J  to  a  governing  crack  propagation  criterion 
in  spite  of  its  underlying  assumptions.  Other  crack  extension  criteria,  closely  associated  with  the  J- 
controlled  growth,  include  crack  opening  angle  (COA),  tearing  modulus  (D,  finite  stretch  and  stability 
index  (A,).  It  is  noted  that  these  extension  criteria  are  all  related  and  permissible  in  treating  a  limited 
amount  of  crack  growth. 

A  leak  before  break  condition  is  expected  for  high  toughness  materials  as  well  as  medium  to  high 
toughness  materials  at  shelf  temperature.  The  large  critical  flaw  sizes  associated  with  these  materials 
suggest  that  extensive  leaking,  as  well  as  general  yielding,  will  precede  attainment  of  a  critical  fracture 
condition.  There  are  simplified  methods  to  compute  mouth  opening  for  pipes  of  these  materials  in 
transition  region;  however,  their  accuracy  remains  to  be  established. 
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2.0  INTRODUCTION 

The  purpose  of  this  work  is  to  study  nuclear  reactor  pipe  failure  due  to  elastic/plastic  fracture. 
The  scope  of  this  study  includes  definition  of  critical  flaw  size-stress  level  conditions  necessary  to 
induce  structural  instability.  Subcritical  growth  mechanisms,  such  as  fatigue  and  stress-corrosion,  as 
well  as  crack  arrest  will  not  be  discussed.  A  primary  objective  here  is  to  assemble  and  to  ultilize  exist¬ 
ing  experimental  data  and  analytical  tools  toward  the  study  of  elastic/plastic  pipe  rupture.  Therefore, 
no  basic  research  in  the  fundamentals  of  elastic/plastic  fracture  mechanics  or  experiments  in  generating 
new  data  are  included.  Only  limited  amount  of  information  is  included  for  the  analytical  methods 
presented  in  this  report.  Details  of  the  formulation  and  justification  should  be  referred  to  the  original 
work  which  are  included  in  the  references.  In  this  report,  vast  amount  of  experimental  data  has  been 
collected  for  later  use;  but  details  of  experimental  techniques  and  variables  are  partially  omitted. 

2.1  LWR  Cracking  Experience 

Cracking  of  LWR  piping  was  first  observed  Dec.  1965,  when  a  leakage  was  found  in  a  6  in.  bypass 
line  of  the  recirculation  loop  in  Dresden  I.  Between  1965  and  1975,  cracks  were  discovered  on  many  4 
in.  diameter  304  s.s.  pipes  for  recirculation  loop  valve  bypass,  and  on  10  in.  diameter  304  s.s.  reactor 
core  spray  lines  (Figs.  1  and  2)  in  six  domestic  Boiling  Water  Reactors.  These  BWRs  are  Dresden  2, 
Quad  City  1  and  2,  Millstation  1,  Peach  Bottom  3  and  Monticello.  Subsequent  investigation  [1]  has 
concluded  these  flaws  were  produced  by  intergranular  stress  corrosion  cracking  (IGSCC)  and  are  due  to 
the  combined  effects  of  stress,  oxidization,  sensitization  and  fatigue  growth.  Subsequent  to  1975, 
IGSCC  has  also  been  found  in  other  location  such  as  reactor-water-clean  up  lines  and  control-rod- 
drive-return  lines.  In  1978,  cracks  in  large  diameter  pipes  (greater  than  20  in  O.D)  were  discovered. 
In  this  case,  extensive  cracking  was  found  on  a  24  in.  diameter  recirculation-inlet-nozzle  safe  end  (Fig. 
3)  at  the  Duane  Arnold  Plant.  Following  intensive  investigation,  NRC’s  pipe  crack  study  group  (2] 
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RECIRCULATION  INLET  NOZZlE 
CARSON  STEEL.  STAINLESS  STEEt  CUB 


Fig.  3.1  — Duane  Arnold  recirculation-inlet-nozzle  safe  end  configuration 


mm  igscc 

•  MEASURED  C«AC«  BERTH 
- ESTIMATED  CRACK  BERTH 

Fig.  3.2— Representation  of  IGSCC  in 
Duane  Arnold  leaking  recircuiation- 
inlet-Nozzle  safe  end 


concluded  that  the  crevice  geometry  at  this  location  and  the  sensitized  material  enhanced  the  IGSCC 
mechanism  even  on  these  larger  pipes. 

The  stress  corrosion  problem  in  PWR  is  not  as  severe  as  that  in  BWR  due  to  low  oxygen  content 
in  the  water  and  fewer  furnace-sensitized  safe  ends.  Therefore,  no  problem  has  been  experienced  in 
PWR  primary  systems.  Nevertheless,  in  the  secondary  system,  inter  and  transgranular  stress  corrosion 
cracking  has  been  observed  in  these  locations  where  relatively  stagnant  boric  acid  solution  are  present. 
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The  impurities  were  introduced  by  safety  injection  and  borated  water  make-up  systems.  At  Arkansas  1 , 
Ginna  and  Surry  1  plants,  the  piping  involved  were  Type  304  s.s.  in  8  in.  and  10  in.  sizes. 

A  detail  complication  of  the  LWR  failure  experience  can  be  found  in  References  3  and  4. 

2.2  Materials,  Geometry  and  Environment 

The  most  commonly  used  materials  in  the  LWR  piping  system  are  Types  304  and  316  austenitic 
stainless  steel  (cast/ wrought).  However,  for  various  reasons,  such  as  intergrrnular  stress  corrosion 
prevention  or  others,  ferrific  steels  such  as  SA-333,  SA-106  and  SA-516  have  also  been  used.  Feed 
water  lines  in  PWR  and  steam  lines  in  BWR  are  typical  examples.  Table  1  is  a  comparison  of  the  ten¬ 
sile  properties  of  these  materials.  In  addition,  low  carbon  stainless  steel  304L  or  306L  and  other  stress 
corrosion  resistant  materials  have  been  recommended  for  piping  applications  [1], 

The  piping  system  in  LWR  is  very  complex  and  a  typical  4-loop  Westinghouse  configuration  is 
shown  in  Fig.  4  [21.  It  can  be  classified  by  its  functional  requirements  or  by  the  material,  geometry  and 
environment  to  which  it  is  subjected.  In  this  report,  since  only  the  consequence  of  the  existence  of  a 
flaw  is  to  be  investigated,  a  pipe  is  referred  by  its  size,  material  composition,  and  the  loading  on  it.  At 
this  point,  it  should  be  noted  that  nearly  all  the  cracks  discovered  to  date  are  located  in  the  weldment 
or  HAZ  where  the  piping  is  connected  to  the  noztle.  These  connections  may  be  to  the  reactor  vessel, 
steam  generator,  feedwater  system  pressurizer,  or  other  components.  Although  the  major  loads  on  the 
system  are  pressure,  thermal  and  mechanical  loads  (e.g.  seismic  and  water  hammering),  the  contribu¬ 
tion  of  the  residual  stress  due  to  the  welding  plays  a  very  important  role  in  initiation  and  propagation 
characteristics  of  cracks.  If  pressure  stress  dictates  the  pipe  failure,  only  axial  cracks  resulting  from 
large  hoop  stress  are  possible.  However,  the  combined  effects  of  material  sensitization,  residual  stress 
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from  welding  processes,  pressure  and  bending  stress  often  initiate  and  propagate  cracks  in  the  circum¬ 
ferential  directions.  In  addition,  because  both  bending  and  residual  stresses  are  self-equilibrating  quan¬ 
tities,  failure  generally  initiates  at  the  surface  where  the  combined  stress  is  the  largest.  In  the  study  of 
pipe  integrity,  both  surface  and  through-wall  cracks  are  of  equal  importance. 

The  operating  conditions  of  a  typical  Pressure  Water  Reactor  (PWR)  is  at  a  pressure  of  2235  psi 
and  at  temperature  of  650°F.  The  Boiling  Water  Reactor  (BWR)  is  operated  at  around  1035  psi  and 
550°F.  The  size  of  these  pipes  are  approximately  30  in.  diameter  for  main  loop  and  steam  generator 
loop,  and  6  to  14  in.  diameter  for  the  other  large  branches. 

2.3  Material  Characteristics  and  Fracture  Toughness 

Toughness  is  one  of  the  material  properties  that  is  essential  to  the  integrity  of  a  structure.  The 
actual  material  toughness  is  both  geometry  dependent  and  temperature  sensitive.  Because  of  the 
microstructural  variation  and  different  constraint  condition  surrounding  a  sharp  flaw,  a  component  can 
have  brittle  failure,  ductile  rupture  or  mixture  of  both. 

To  illustrate  this  phenomeon,  a  schematic  Dynamic  Tearing  (DT)  test  result  of  fixed  thickness 
specimens  are  shown  in  Fig.  5.  At  the  Nil  Ductility  Transition  (NDT)  temperature,  the  fracture  is  brit¬ 
tle  and  shows  a  flat,  featureless  surface.  A  rapid  increase  in  fracture  energy  is  recorded  at  temperatures 
above  NDT  as  more  ductility  is  developed  which  is  evidenced  by  increase  in  lateral  contraction  and 
development  of  shear  lips.  As  the  temperature  exceeds  shelf  temperature,  brittle  cleavage  appearance 
is  replaced  by  ductile  dimple  type  failure  and  there  is  no  further  increase  in  fracture  resistance  above 
this  point.  The  basics  of  this  fracture  energy/temperature  relationship  can  be  explained  from  the 
differences  in  microstructural  failure  modes.  Brittle  failure  at  temperatures  considerably  below  the  shelf 
involves  pure  cleavage  of  the  individual  grains  and  is  a  high-speed  process.  However,  upon  increasing 
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DT 

ENERGY 


Fig.  5— DT  test  transitions  to  various  levels  of  shelf  fracture  toughness.  Note  that  with  a  decrease  in  shelf  level  frac¬ 
ture  energy  there  is  a  corresponding  change  from  fractures  with  large  lateral  contraction  to  flat  fractures  with  nil  con¬ 
traction  features.  The  decrease  in  shelf  energy  marks  a  transition  from  plastic  (plane  stress)  to  elastic  (plane  strain) 
fracture  conditions. 


in  temperature,  cleavage  separation  of  individual  grains  competes  with  slip  processes.  More  energy  is 
requried  for  attaining  higher  stress  needed  for  cleavage  because  more  strain  is  required  for  elevating 
flow  stress  to  the  level  of  the  cleavage  stress.  Macroscopically,  in  the  transition  region,  the  increase  in 
fracture  toughness  is  accompanied  by  the  development  of  shear  lips  at  the  fracture  surface.  Finally,  at 
the  upper  shelf  temperature,  fracture  process  is  defined  by  microvoid  coalescence;  when  small  voids 
between  grains,  or  of  inclusions  or  impurities,  are  opened  and  the  metal  bridges  between  these  voids 
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are  stretched  as  tiny  tensile  specimens  which  finally  rupture  in  a  progressive  ductile  mode.  The  transi¬ 
tion  temperature  and  the  sharpness  of  transition  region  depends  on  the  specimen  thickness  (Fig.  6). 
Generally,  the  thinner  specimen  exhibits  lower  transition  temperatures. 


Temperature ,  F 


Fig.  6— Effect  of  plate  thickness  on  DWTT  and  Charpy  shear  area  results 


The  validity  of  linear  elastic  fracture  mechanics  (LEFM)  is  limited  to  small  scale  yielding.  When 
the  size  of  plastic  zone  developed  in  the  specimen  exceeds  a  prescribed  limit  (Fig.  7)  due  to  changes  in 
service  loading  or  increase  in  temperature,  Kk  (6l  or  K,d  looses  its  usefulness.  The  region  of  valid 
LEFM  toughness  determination  is  restricted  to  the  temperature  (T)  and  energy  level  (S)  limits  indi¬ 
cated  in  Fig.  8.  Fortunately,  the  J-integral  [7-16]  approach  has  extended  the  material  toughness  meas¬ 
urement  to  upper  shelf  region.  The  development  of  J-integral  concept  and  testing  technique  not  only 
has  reduced  K  testing  effort  by  using  smaller  specimens  but  also  has  enabled  structural  engineers  to 
characterize  a  flawed  body  subjected  to  large  scale  yielding.  However,  the  validity  of  J-integral  approach 
is  limited  theoretically  due  to  its  basic  assumptions  [7]  such  as  infinitesimal  deformation,  deformation 
theory  of  plasticity  etc. 
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Fig.  7— Relationships  of  elastic  and  plastic  stress  fields  to  the  plastic  zone  at 
crack  lips  for  the  case  of  plane  strain  constraint.  As  plastic  relaxation  is 
developed  (large  plastic  zone  and  crack  tip  blunting),  the  elastic  stress  fields  are 
replaced  by  plastic  strain  fields.  Elastic  stress  field  K  definitions  are  not  possible 
for  these  conditions. 


Under  normal  reactor  operating  conditions,  materials  that  have  been  used  for  nuclear  piping  con¬ 
struction  can  be  identified  as  those  exhibiting  no  transition  temperature  and  those  having  a  definite 
transition  temperature.  Austenitic  stainless  steel,  304  and  316  s.s.,  are  examples  of  the  first  kind  and 
ferritic  steels  are  the  second  type.  Therefore,  austenitic  stainless  steels  are  materials  having  very  high 
toughness  and  the  structural  failure  is  generally  related  to  limit  load  conditions.  Figure  9  illustrates  the 
J-R  curves  of  304  and  316  s.s.  at  room  temperature  as  well  as  600°F  [17].  The  J,c  (critical  initiation 
value)  ranges  from  5000  in.  lb/in.2  at  RT  to  3500  in.  lb/in.2  at  600°F.  Using  the  Klc-Jlc  relationship: 


*/?- 


EJi, 


\-v2 


where  E  is  the  modulus  and  v  is  Poisson's  ratio,  the  critical  stress  intensity  factors  K/c  are  375  ksi  Vim 
and  300  ksi  Vim,  respectively.  With  such  high  toughness  values,  it  is  obvious  that  these  materials  have 
high  tolerance  against  brittle  fracture,  and  the  structural  failure  must  be  controlled  by  ductile  tearing. 


12 


NRL  MEMORANDUM  REPORT  4259 


When  attention  is  directed  to  the  materials  of  second  type  which  have  lower  toughness  and  exhi¬ 
bit  a  definite  transition  temperature,  the  structure  failure  is  no  longer  controlled  by  ductile  tearing 
alone.  Depending  on  the  operating  temperature  regime  of  the  material  and  the  geometry  of  the  struc¬ 
ture  (thickness,  pipe  size  and  surface  flaw  or  through-wall  flaw),  brittle  fracture  and/or  ductile  tearing  is 
possible.  Generally,  the  material  toughness  characterization  tests  and  structural  failure  analyses  of 
these  types  of  materials  are  more  complicated  and  less  conclusive. 

3.0  COMPUTATIONAL  METHODS  OF  NUCLEAR  REACTOR  PIPE  RUPTURE 

Typically,  the  pipe  rupture  study  can  be  classified  as  near  field  or  far  field  study.  The  terms,  near 
and  far,  denote  the  relative  size  of  the  region  around  the  crack  tip  where  the  analysis  is  employed.  In 
the  far  field  computational  method,  a  considerable  amount  of  effort  is  directed  to  model  the  overall 
structural  response  of  piping  system,  and  the  crack  tip  behavior  enters  into  the  analysis  only  as  a  simple 
fracture  criterion.  On  the  other  hand,  near  field  study  deals  with  the  materials  and  structural  response 
locally  near  the  cracked  region.  The  overall  structural  response  affects  the  local  stability  only  in  the 
compliance  formulation. 

In  the  far  field  computation,  the  whole  piping  system  is  generally  modeled  by  3-D  solid  or  shell 
finite  elements.  The  dynamic  behavior,  produced  by  an  advancing  crack,  of  the  piping  system  and  the 
fluid  inside  are  considered.  Because  of  the  computational  complexity  in  dealing  with  the  overall  struc¬ 
tural  dynamics  of  the  piping  system,  only  the  simple  crack  extension  criterion  can  be  included  in  the 
iterative  process  without  prohibitive  computational  time.  The  commonly  used  criteria  are  either  max¬ 
imum  strain  or  K/c. 

In  the  LEFM  analysis,  toughness  can  be  considered  as  a  material  property.  The  stability  of  a 
flawed  structure  is  conservatively  assured  when  the  applied  stress  intensity  factor  is  less  than  the  critical 
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value,  Kk.  This  is  the  basic  philosophy  that  the  ASME  boiler  and  pressure  vessel  code  has  adopted. 
Because  the  structural  stability  is  measured  by  a  near  field  parameter  A,  this  type  of  analysis  can  be 
considered  as  a  near  field  analysis.  Another  example  of  the  near  field  approach  is  the  fracture  mode 
transition  temperature  method,  where  the  operating  temperature  is  used  to  control  structural  failure 
against  fracture  initiation  and  propagation  of  existing  flaws.  However,  when  cracks  advance  into 
material  resulting  in  large  deformation,  stability  against  fracture  may  depend  on  material  as  well  as 
geometry  and  loading  conditions.  In  this  case,  better  failure  description  is  needed.  The  J-integral 
approach  is  rigorous  in  defining  fracture  initiation  under  gross  plastic  yielding.  Evidence  to  date  indi¬ 
cates  that  the  J-R  resistance  curve  may  be  useful  in  studying  crack  propagation  under  monotonic  load¬ 
ing.  If  one  accepts  this  assumption,  the  J-R  resistance  curve  approach  can  also  be  considered  as  a  near 
field  analysis  method. 

Both  near  field  and  far  field  methods  assume  that  the  local  instability  criteria  are  dependent  only 
on  the  material.  However,  the  structural  geometry  and  the  loading  are  needed  in  the 
structural/ material  response  computation.  To  make  this  complicated  problem  more  trackable,  much 
effort  has  been  directed  to  develop  a  simple  analytical  formulation  which  considers  the  material  proper¬ 
ties  as  well  as  structural  geometry  and  loading  conditions.  The  technique  used  in  generating  these  for¬ 
mulations  is  generally  semi-empirical  in  nature.  Starting  from  a  known  solution  for  a  flawed  plate, 
modifications  to  account  for  plasticity  and  geometry  effects  are  first  incorporated.  The  final  simplified 
solution  is  derived  by  verification  and  adjustment  of  the  modified  formulation  using  a  vast  amount  of 
experimental  data. 

In  the  next  several  sections,  detailed  discussion  of  each  computational  method  will  be  presented. 
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3.1  Gross  Structural  Response 

This  type  of  computational  method  is  far  field  study  because  the  analysis  includes  a  very  large 
region  of  the  piping  system.  In  addition,  the  computational  emphasis  is  on  the  overall  response  of  the 
structural  rather  than  local  region  surrounding  an  existing  crack.  With  improved  numerical  computa¬ 
tion  capabilities,  pipe  rupture  can  be  very  easily  modeled  by  finite  element  or  finite  difference  computer 
programs.  In  these  analysis,  typical  shell  elements  are  used  to  model  the  piping.  Elastic,  elastic/plastic, 
or  even  viscoelastic  constitutive  behavior  can  be  employed.  In  addition,  the  dynamics  of  crack  opening, 
coupled  with  the  escape  of  the  internal  fluid,  can  also  be  included.  Depending  on  the  requirement,  a 
very  costly  program  may  be  developed  to  model  some  postulated  event.  One  of  the  key  ingredients  in 
simulating  the  piping  rupture  is  the  requirement  of  an  adequate  fracture  criterion  such  as  maximum 
strain  or  Kk.  References  18  &  22  are  typical  research  results  of  this  kind.  However,  the  accuracy  of 
these  analytical  results  depends  heavily  on  a  poorly-defined  quantity,  viz.,  the  fracture  toughness  of  the 
material.  When  brittle  failure  is  observed,  Klc  or  K]d  may  be  adequate.  For  a  pipe  that  fails  plastically 
or  in  a  mixed  brittle/ductile  failure  mode,  the  analytical  prediction  can  be  misleading.  In  addition,  the 
numerical  results  of  one  costly  numerical  analysis  may  not  be  extrapolated  to  other  cases  where  the 
prohibitive  costs  prevent  additional  runs.  The  need  for  this  type  of  analysis  is  apparent  when  dynamic 
fracture  (static  initiation/dynamic  propagation)  governs  the  fracture  processes.  However,  based  on  the 
documented  nuclear  reactor  piping  field  failure  experience,  there  is  very  little  evidence  that  a  pipe 
failure  due  to  dynamic  propagation  is  likely.  The  primary  reason  for  a  non-propagating  crack  is  the 
fluid  (water)  does  not  enhance  dynamic  growth  in  piping.  In  general,  gross  structureal  analysis  is  not  a 
very  widely  accepted  method  in  pipe  rupture  analysis. 

3.2  Semi— Empirical  Methods 

Before  an  acceptable  inelastic  pipe  rupture  analytical  method  is  developed,  the  most  logical  means 
to  study  the  problem  is  to  rely  on  known  technology.  In  this  case,  the  existing  technology  is  linear 
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elastic  fracture  mechanics.  However,  one  immediately  faced  with  the  task  of  circumventing  the 
inherent  limitations  built  within  LEFM.  The  test  specimen  thickness  requirement  is  one  example.  The 
constraint  requirement  imposed  by  ASTM  E-399  with  respect  to  determination  of  the  plane  strain  criti¬ 
cal  stress  intensity  factor,  Klc,  is 

B  >  2.5  \^~  .  (1) 

I*,* 

For  materials  having  yield  strength  <rys ,  the  thickness  (B)  must  be  in  excess  of  the  actual  reactor  piping 
thickness  to  have  an  adequate  K,c  measurement.  Since  K,c  is  very  sensitive  to  thickness  variation, 
experimentally  generated  K,c  values  can  be  misleading.  Another  important  consideration  in  adopting 
the  LEFM  approach  is  extensive  plastic  yielding  during  the  failure  process.  If  pipe  failure  is  due  to 
gross  plastic  yielding,  the  elastically-generated  fracture  toughness  value  becomes  meaningless. 

To  account  for  the  LEFM  limitations  as  well  as  the  geometry  and  loading  conditions  on  a  pipe, 
semi-empirical  methods  have  been  proposed.  In  this  respect,  either  a  pseudo-toughness  value,  Kc,  or 
flow  stress  theory  is  applicable,  depending  upon  the  fracture  mode  (brittle,  brittle/ ductile,  ductile) .  The 
pseudo-toughness  value,  Kc,  is  an  equivalent  critical  toughness  parameter  for  a  piping  application.  It  is 
derived  by  modifying  elastic  fracture  mechanics  solutions  of  a  flawed  plate  through  curvature 
correction,  plasticity  consideration,  and  extensive  experimental  pipe  rupture  data  correlation.  When 
pipe  failure  is  dictated  by  gross  plastic  yielding,  a  flow  stress  theory  is  derived  by  assuming  that  pipe 
instability  is  governed  by  the  limit  load  of  the  pipe. 

It  will  become  apparent  in  subsequent  sections  that  these  semi-empirical  methods  are  easy  to  use 
but  do  not  assure  good  correlation  with  the  actual  pipe  test  results. 
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3.2.1  Axial  Cracks 

For  an  axial  crack  in  gas  transmission  line  pipe,  equation  (2)  has  been  proposed  [23]  to  compute 


Kc. 


,  rrcrrl  c 2 

K ? - -f  1  +  1.61  -s- 

f  cos  0  Rt 


4 -K 


K} 


*c<T(i 
cos  0 


(M1) 


4  -  K 


where 


(2) 

(2a) 


Kc  —  critical  stress  intensity,  ksix/TrT 
c  *=  half  axial  through-wall  crack  length,  in. 
R  —  average  radius  of  vessel,  in. 
t  —  wall  thickness,  in. 


<r  h  -  — ,  nominal  hoop  stress  at  failure,  ksi 

P  —  vessel  failure  pressure,  ksig 

r  —  inside  radius  of  vessel,  in. 


<rc  —  failure  stress  for  unflawed  vessel,  ksi 
(M<rh  was  used  for  this  value) 

K  —  (3  -  4v)  plane  strain;  (3  -  v)/(l  +  v)  plane  stress 


19 


CHANG,  NAKAGAKI,  GRIFFIS  AND  MASUMURA 


v  —  Poisson's  ratio 


M  —  stress  magnification  factor  for  an  axial  through  crack  in  a  cylinder  which  is  a  function  of  A 
(see  Figure  10;  A/exacl  exact  was  used  herein) 


For  axial  cracks  in  an  intermediate  wall  thickness  vessel  or  pipe,  Folias  [241  and  Goodier  &  Field 
1251  proposed  equation  (3)  for  low  to  medium  toughness  materials  with  relatively  long  crack.  Hahn 
[261  derived  equation  (4)  for  high  toughness  materials  with  short  cracks. 


a*  is  the  flow  stress  of  the  material,  which  according  to  Hahn  can  be  taken  as  1.04  o-,  +  10.0 
(ksi)  or  with  less  accuracy  by  0.51  (<ry  +  <rM) 

cr^,  —  yield  stress,  ksi 


<tu  —  tensile  stress,  ksi 

f  ,2  "2 

M  —  II  +  1.61  —  or  from  Figure  10  for  a  better  estimate 

tr  *  -  trhM  (4) 

[  c2  1/2 

where  M  —  1  +  1.61  ~  or  A/„acl  from  Figure  10  for  a  better  estimate.  The  relationship  between 

equations  (3)  and  (4)  is  illustrated  in  Fig.  11.  It  is  seen  that  equation  (4)  is  the  upper  limit  for  the 
high  toughness  (ductile)  behavior. 
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10  12 


2c  A/Rt 


Fig.  1 1 —Dimensionless  failure  curves  for  Eq.  (3)  and  (4) 


To  describe  the  failure  of  low  to  medium  toughness  materials,  Newman  [27]  proposed  a  two- 
parameter  fracture  criterion  to  include  material  yielding  as  well  as  brittle  failure.  These  two  parameters, 
m  and  Kf,  are  computed  from  laboratory  test  results  (Attests)  and  are  related  to  the  net  section  failure 
stress  S„,  of  a  specimen  having  half  crack  length  c,  width  W,  and  ultimate  strength  <ru. 

N  S  H  N  N  S 

— 

1-1  (-1  (-1  ,-|  O’u 

m  V  s  )2  a  N  Ts  r 

I  —  I H-  -==- 

i-\  lau  j-1  /-I  <^u 


N 

I  Kl 

Kf - — - 

I  Ku-m  ~ 

i— 1  i—  I  t'Tu 

where  Kte  is  the  computed  elastic  stress  intensity  based  on  gross  section  stress  5  to  failure  equation 


sec  (irc/w). 
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For  plane  strain  behavior  m  approaches  zero  Kf  — •  Kle  —  Klc  (the  plane  strain  fracture  toughness). 
For  notch  insensitive  materials,  m  becomes  unitv  so  Kf  relates  failure  to  the  ultimate  tensile  strength. 
In  the  range  covered  by  mixed  mode  and  plane  stress  fracture,  the  failure  stress  is  a  function  of  both 
Kf  and  m.  To  account  for  the  curvature  effect  of  a  cylindrical  structure,  Adams  [28]  proposed  a  revised 
failure  equation  which  gives  the  failure  stress,  5cal,  as 


Seal  ~ 


C,-Jnc  —  sec  ( nc/w )  + 


mK 


f 


As 


An 


(8) 


where  Ag  and  A„  are  cross  sectional  area  and  net  section  area,  respectively.  A  curvature  correction  fac¬ 
tor,  Cc,  is  defined  by: 


Cc  =  0.614  +  0.48A  +  0.386e(_l  25*\ 

where 


(9) 


A  =  [12(1 -*2)1/4c]/V7?7.  (10) 

Figure  12  [29]  illustrates  the  comparison  of  equations  (2),  (3)  and  (8)  for  a  42  in.,  3  in.  wall  thickness 
pipe.  Approximations  made  in  equation  (8)  were  m  ~  1.0  and  nc/w  0.0.  It  is  noted  that  for  Kc  less 
than  200  ksi  Vim,  all  three  equations  have  similar  toughness  ( Kc )  and  flaw  size  (2c)  relationship.  This 
observations  is  encouraging  because  most  reactor  piping  materials,  displaying  brittle  to  ductile  transition 
iperature  response,  exhibit  toughness  ranges  from  50  to  200  KSI  VTnT 


Other  semi-empirical  equations  to  define  piping  toughness  have  been  described  by  Folias  [31]  and 
Quirk  [30],  However,  a  detailed  assessment  program  using  A  106  B  pipe  experiments  [31]  indicated 
that  equations  (3)  and  (4)  are  the  most  suitable  equations  to  be  used.  There  is  evidence,  however,  that 
the  scatter  of  the  Kc  computation  may  be  much  larger  than  that  indicated  in  Ref.  [31]  when  materials 
or  different  pipe  sizes  are  considered.  It  is  concluded  that  the  pseudo-fracture  toughness,  Kc,  can  be 
used  as  a  toughness  indication  and  it  is  accurate  under  specific  conditions.  More  discussions  on  this 
point  will  be  followed  in  a  subsequent  portion  of  the  text. 
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Equation  (8). Ref.  28 
9  \ 


□ 


25  A 


*  20  I 


- □ 

A 


Equation  (2),  Kef. 23 


B  O 
o 


-p 

•H 

O 


o 

C\J 


15-1 


10 


8 


Equation  31,  Kef. 25 


NRL  MEMORANDUM  REPORT  4259 


3.2.2.  Circumferential  Cracks 


Circumferential  cracks  can  only  be  initiated  and  propagated  by  large  axial  stress.  Therefore,  the 
necessary  condition  associated  with  these  types  of  cracks  is  the  presence  of  high  secondary  stresses 
(thermal,  bending,  expansion  and  residual  stresses).  Although  most  pipe  failures  in  the  field  are  due  to 
circumferential  flaws  (5),  no  rigorous  analytical  method  is  available  to  date.  Laboratory  tests  on 
medium  to  high  toughness  materials  [32,33,34,35]  revealed  that  for  short  circumferential  flaws  the  pro¬ 
pagation  was  axial  from  the  crack  tips.  This  result  indicated  that  prototype  tests  were  not  good  simula¬ 
tions  to  account  for  the  complicated  stress  state  and  microscopic  material  degradation  of  the  materials  in 
service.  However,  when  longer  circumferential  flaws  were  studied,  circumferential  propagation  was 
observed.  In  this  case,  tests  on  high  toughness  materials  seem  to  correlate  well  with  an  ultimate 
strength  theory.  This  theory  suggests  that  pipe  rupture  is  produced  when  the  axial  stress,  due  to  pres¬ 
sure  and  bending,  reaches  the  material  ultimate  strength.  The  bending  stress  in  the  pipe  is  derived  by 
accounting  for  the  shifting  of  the  section’s  neutral  axis  due  to  the  presence  of  a  flaw.  One  of  the  for¬ 
mulations  [3]  gives  the  corresponding  limit  moment  (Fig.  13)  as: 


where 


M  = 


4 (it  —  a)2R2t2<rl  -  n2R2P2 
2(tt  -  a)2R2ta-0 


[/?o  (2 cos/3  -  sina)] 


(11) 


<t0  —  flow  stress 

Rm  =  mean  radius 

Ri  —  inside  radius 

R o  “  outside  radius 

T  “  thickness 

P  -  internal  pressure 
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Another  method  that  may  be  useful  is  the  solution  developed  by  Erdogan  and  Kibier  [36).  In 
their  work,  stress  intensity  factors  have  been  computed  as  a  function  of  the  parameter  A  (Eq.  10).  The 
cylinder  is  subjected  to  both  axial  and  bending  loads.  However,  lack  of  a  plastic  zone  correction  factor 
renders  their  formulation  somewhat  unrealistic. 

3.3  Semi-Empirical  Methods  Correlation 

The  adequacy  of  the  proposed  semi-empirical  methods  are  highly  dependent  on  the  material  being 
considered.  Namely,  the  pseudo-toughness  method,  Kc ,  is  for  medium  to  high  toughness  materials, 
whereas  the  flow  stress  theory  is  suited  for  high  toughness  material  applications.  In  the  following  dis¬ 
cussion,  correlations  between  analysis  methods  and  experimental  test  data  are  given  for  each  type  of 
material.  Both  axial  and  circumferential  cracks  are  addressed.  It  will  become  apparent  that  both 
methods  can  not  be  used  indiscriminately.  Lack  of  good  correlation  between  test  data  and  projections 
according  to  the  pseudo-toughness,  Kc,  method  raises  doubts  about  the  adequacy  of  this  approach.  On 
the  other  hand,  even  with  limited  success,  the  flow  stress  theory  is  still  a  questionable  method  to 
predict  inelastic  pipe  rupture.  Specifically,  because  the  flow  stress  theory  is  toughness  independent  and 
most  of  all,  geometry  independent.  It  is  felt  that  the  flow  stress  theory  can  be  used  with  success  for 
most  of  the  cases.  But  before  one  adapts  its  methodology,  careful  consideration  should  be  given  to  the 
specific  piping  system  studied  and  the  assumptions  and  limitations  of  the  theory. 

3.3.1  High  Toughness  Materials 

High  toughness  materials,  such  as  304  and  316  stainless  steels  exhibit  no  brittle  to  ductile  transi¬ 
tion  temperature  and  generally  fail  by  ductile  rupture.  The  measured  high  toughness  (Fig.  9)  seems  to 
rule  out  the  possibility  of  brittle  fracture  initiation.  However,  cracks  have  been  found  in  these  materi¬ 
als  due  to  various  causes.  When  this  material  is  subjected  to  heat  treatment  between  800°F  and 
1200°F,  chromium  is  depleted  from  the  matrix  by  precipitation  at  the  grain  boundaries  in  the  from  of 
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chromium  carbide.  The  material,  having  undergone  this  metallurgical  change,  is  said  to  be  sensitized. 
When  the  proper  agent,  such  as  oxidizing  element,  is  introduced  to  the  sensitized  material,  stress  corro¬ 
sion  cracking  can  be  developed  under  imposed  thermal  and  mechanical  loads.  In  this  case,  the  objec¬ 
tive  is  not  only  to  devise  techniques  to  reduce  stress  corrosion  possibility,  but  also  to  evaluate  the  pip¬ 
ing  stability  due  to  the  presence  of  critical  flaws  in  various  orientations. 

It  has  been  suggested  [37,38,39,40]  that  the  flow  stress  theory  (Eq.  4)  may  be  an  adequate 
method  in  predicting  ductile  fracture  failure  of  high  toughness  materials.  As  suggested  in  the  following 
sections,  the  limited  experimental  results  appear  to  correlate  well  with  the  theory.  However,  because 
no  crack-tip  parameter  is  considered  in  the  formulation,  it  is  premature  to  conclude  that  this  theory  is 
applicable  regardless  of  piping  geometry  or  applied  loads.  Also,  Tada  and  Paris  [65]  have  demonstrated, 
using  J  integral  techniques,  the  importance  of  pipe  length/diameter  (L/R)  ratio  in  a  ductile  stability 
analysis  of  a  circumferential  crack  in  reactor  piping. 

3.3. 1.1  Axial  Cracks 


Battelle  [37,38]  carried  out  4  experiments  on  24  in.  diameter  316  stainless  steel  pipes  of  1.5  in. 
thickness.  The  length  of  these  pipes  was  between  8  to  22  ft.  The  pipes  were  heated  between  470- 
680°F  and  pressurized  to  5000  psi  to  induce  failure.  Both  surface  cracks  and  through-wall  cracks  were 
introduced  in  the  specimens  (Fig.  14).  Using  equation  (4),  corrected  for  reductions  in  area  due  to  exis- 
tance  of  surface  flaw,  the  flow  stress  equation  becomes  [37] 


t!d  -  1 


(12) 


*  ”  t/d-  1/M 

where  <rh  is  the  hoop  stress  at  failure;  <r*  is  the  flow  stress;  M  is  the  stress  magnification  factor  as 
defined  in  Fig.  10;  and,  rand  d  are  pipe  thickness  and  flaw  depth,  respectively.  When  very  long  flaws 
are  considered,  l/M  approaches  zero  and  the  equation  (12)  is  reduced  to 


<rh  -  cr*(t-d)/l. 


(13) 
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Electric  Fracture 


In  this  case  failure  stress  is  proportional  to  the  remaining  ligament  thickness.  On  the  other  hand,  when 
deep  flaws  are  considered,  equations  (12)  reduces  to  the  flow  stress  equation  (4)  for  through-wall 
cracks.  Figure  15  illustrates  the  correlation  between  test  points  and  computed  values;  and  it  is  evident 
that  good  agreement  was  obtained.  Since  316  s.s  is  a  very  high  toughness  material,  even  at  room  tem¬ 
perature,  gross  plastic  deformation  before  failure  is  expected.  Therefore,  for  the  pipe  geometry  stu¬ 
died,  the  flow  stress  criterion  is  adequate. 

3.3.1.2  Circumferential  Cracks 

To  investigate  the  austentic  stainless  steel  toughness  against  circumferential  flaws,  Battelle  [39,40] 
undertook  a  series  of  experiments  performed  on  type  304  stainless  steel  plate  and  pipes.  Flat  plate 
specimens  containing  center  cracks  were  used  to  evaluate  the  effect  of  HAZ  (sensitization)  and  crack 
tip  sharpness  on  the  gross  behavior  of  the  material.  It  was  found  that  extensive  blunting  practically 
overshadows  the  effect  of  initial  crack  tip  geometry  on  the  final  instability.  Also,  the  propagation 
characteristics  were  very  similar  in  HAZ  and  base  metal.  Because  the  exhibited  high  toughness  values 
of  304  s.s.,  the  results  from  the  plate  tests  suggested  that  the  flow  stress  theory  is  an  adequate  criterion 


Fig.  15— Surface  flaw  behavior  in  24  x  1  50  Type  316  stainless  steel  pipe 
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to  predict  collapse  loads.  Ultilizing  the  computed  flow  stress  values  from  the  plate  tests  and  equation 
(11),  the  limiting  moment  versus  flaw  length  relationship  is  plotted  for  internal  pressures  of  1050  psi 
and  2500  psi  (Fig.  16).  Experimental  points  were  obtained  from  full-scale  pipe  experiments  on  two  4- 
in. -diameter  schedule  80  type  304  stainless  steel  pipes.  Initial  circumferential  flaws  subtended  arcs  over 
135°  and  75.8°,  and  testing  was  performed  at  a  temperature  of  approximately  3°C.  A  good  correlation 
between  test  results  and  flow  stress  projections  was  observed. 

For  the  case  of  surface  cracks,  equation  (12)  can  be  modified  by  including  the  ligament  area  in 
the  derivation  of  the  cross  sectional  area  characteristics  [38]. 

3.3.2  Medium  to  High  Toughness  Materials 

In  nuclear  piping  applications,  this  type  of  material  exhibits  a  brittle  to  ductile  transition  tempera¬ 
ture.  It  can  be  considered  as  high  toughness  material  when  it  is  operated  at  or  above  the  shelf  tempera¬ 
ture.  Therefore,  the  failure  modes,  and  consequently  the  analysis  methods,  are  temperature  dependent. 
For  this  reason,  more  sub-scale  and  full-scale  experiments  have  been  performed  on  piping  and  vessels 
of  this  type  of  material.  It  is  evident  from  the  discussion  presented  in  following  sections,  that  there  are 
a  number  of  technical  questions  which  remain  unanswered.  The  limited  evidence  to  date  indicates  that 
when  these  materials  are  operating  at  or  above  shelf  temperature  flow  stress  theory  (Eq.  4)  is  adequate 
in  collapse  load  prediction.  Nevertheless,  the  inherent  simplicity  in  the  flow  stress  theory  (Section 
3.3.1)  prevents  one  from  adopting  this  theory  for  arbitrary  pipe  geometry  and  the  imposed  system  load¬ 
ings.  When  the  operating  temperature  drops  into  transition  region,  the  pseudo-fracture  toughness,  Kc, 
method  may  be  used  for  selected  axial  flaws;  however  no  adequate  analytical  solution  currently  exists 
for  circumferential  cracks.  Furthermore,  the  poor  correlation  between  Kc  and  test  data  raises  the  ques¬ 
tion  of  the  general  adequacy  of  the  pseudo-toughness  theory  for  piping  analysis.  An  additional  practical 
consideration  that  contributes  towards  the  difficulty  in  analyzing  these  materials  is  the  scatter  in  material 
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3.3.2. 1  Axial  Cracks 

For  materials  containing  axial  cracks  and  operating  at  shelf  temperatures,  the  flow  stress  theory 
(Eq.  4)  seems  to  be  appropriate  for  the  selected  cases  investigated.  Battelle’s  [37,38]  work  on  A106  B 
carbon  steel  is  a  typical  example  of  on-the-shelf  behavior.  A  total  22  experiments  were  conducted  on 
12.75  in.  diameter  and  24  in.  diameter  pipes  containing  surface  and  through-wall  cracks.  These  pipes 
were  heated  above  450°F  under  internal  pressure  (Table  2).  According  to  the  Charpy-V-Notch  test 
data  contained  in  the  Appendix  A  [24],  450°F  may  be  considered  to  be  a  shelf-level  temperature.  Fig¬ 
ure  17  illustrates  the  excellent  agreement  between  test  point  and  the  flow  stress  criterion,  (Eq.  4). 

The  on-the-shelf  behavior  of  surface  flaws  on  A 106  B  steel  pipes  can  also  be  assessed  from  other 
experimental  programs.  Table  3  contains  the  test  conditions  of  BatlelJe’s  work  on  surface  flaws  [37], 
To  investigate  the  validity  of  flow  stress  criterion,  failure  stress  to  flow  stress  ratios  (o-A/o-*)  are  com¬ 
puted  for  experimental  results.  Corresponding  theoretical  <r  Jo  *  values  are  also  computed  from  Eq.  12 
for  the  d/t  rations  tested.  In  Fig.  18  a  perfect  correlation  line  at  45  degrees  is  shown,  and  the  excellent 
agreement  between  theory  and  experiment  is  noted. 

For  ferritic  piping  materials  operating  in  the  brittle  to  ductile  transition  region,  experimental  data 
have  been  generated  to  characterize  the  fracture  response.  General  trends  have  been  observed  for 
specific  materials  and  specific  pipe  geometries.  Figure  19  illustrates  the  failure  stress  vs.  temperature 
relationships  generated  from  UKAEA  data  [5,41]  (Tables  4,  5)  on  0.36%  carbon  steel.  It  is  noted  that 
6-in. -flaw  failure  stress  curve  crosses  the  yield  strength  curve  at  135°F,  whereas  larger  flaws  push  the 
crossover  point  to  higher  temperatures.  Because  smaller  flaws  produce  full  ductile  rupture  at  lower 
temperature,  this  implies  that  smaller  flawed  structures  have  lower  transition  temperatures.  GE  results 
[5,42]  on  schedule  40  A 106  B  pipe  at  room  temperature  (Table  6)  indicated  that  larger  pipe  have  higher 


TABLE  2  CRITICAL  CRACK  EXPERIMENTS  ON  CARSON  STEEL  PIPE  WITH  THROUGH-WALL  FLAWS 


CIO  21  605  2.5  ?9.0  *2.3  74.0  6.375  0.710  1.  10  166  42.9  12.0 

Averages  42.8  74.0 


SOLID  POINTS  ARE  EXPERIMENTAL  FAILURE 
POINTS.  ALL  FRACTURES  IN  THE  EXPERIMENTS 
WERE  COMPLETELY  DUCTILE. 


Comparison  of  experimental  critical  flaw  si/e  data  with  calculated  behavior  for  three  pipe  sizes 
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(b) (ou+  a  y)  /2.4. 

(c)  Predicted  stresses  are  based  or.  a*  (com  tensile  test  results. 

(d)  This  leaked  because  G$at-  was  too  low. 
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TABLE  4.  VKAEA  AXIAL  UMOUCH'UAU.  FLAW  MIA  OH  JTTtL  CTUHOlkCAL  VUttU 


i 


TABLE  6.  THROUGH-WALL  FLAW  DATA  OK  A10CB  PIPE  FROM  REYNOLDS 


(c)  Oc  was  ta/en  as  (CTU  ♦  oy)/2.4 
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tolerance  to  the  existence  of  larger  axial  flaws  (Fig.  20).  To  date,  macroscopic  observations  and  analyti¬ 
cal  solutions  have  been  able  to  establish  flawed  structural  behavior  with  limited  success.  Specifically,  a 
non-dimensional  parameter  0  seems  to  be  one  of  the  variables  that  can  be  ultilized  to  correlate  the 
laboratory  results.  The  shape  parameter,  0,  is  defined  by: 

*-7k 

where 

2c  —  flaw  length 

R  “  pipe  radius 

t  =  pipe  wall  thickness 

To  test  this  assumption,  a  correlation  between  failure  stress  and  0  is  made  on  Fig.  21  [371.  On 

this  graph  (all  test  data  are  reproduced  in  Appendix  A),  the  flow  stress  theory  (Eq.  4)  is  shown  as  solid 
lines  and  McDermott’s  [43]  limit  analysis  is  shown  as  dotted  lines.  A  general  trend  is  seen,  but  the 
scattered  data  points  signal  the  importance  of  other  variables  which  have  not  been  incorporated  into  the 
definition  of  0. 

As  stated  previously,  the  uncertainty  in  as-received  material  properties,  flaw  geometry,  pipe  size 
and  pipe  geometry  play  a  very  important  role  in  the  precise  location  of  the  brittle  to  ductile  transition 
temperature  of  a  specific  structure.  When  the  temperature  is  below  NDT,  LEFM  is  adequate  and  when 
the  shelf  temperature  is  reached,  flow  stress  criterion  may  be  appropriate.  However,  in  the  transition 
region,  mixed  mode  of  failure  occurs.  The  degree  of  ductility  associated  with  failure  is  affected  by  vari¬ 
ous  factors  and  it  is  very  difficult  to  express  fracture  response  in  the  framework  of  present 
elastic/ plastic  fracture  mechanics  technology. 

Among  the  existing  criteria,  pseudo-toughness  ( Kc )  and  flow  stress  (<r*)  theories  seem  to  be 
favored  by  researchers  because  of  their  simplicity.  However,  if  a  material  fails  in  mixed  mode  manner. 


Failure  Stress,  ksi 
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application  of  the  flow  stress  theory  is  tenuous  due  to  its  underlying  assumption  of  full  ductility.  On 
the  other  hand,  the  pseudo-toughness  formulation  has  the  capability  to  rigorously  treat  brittle  failure 
and  also  incorporates  a  plasticity  correction.  Unfortunately,  the  semi-empirical  nature  of  K ,  requires 
extensive  experimental  data  to  justify  its  validity  and  also  limits  its  ability  to  extrapolate  to  different 
ma'erials,  structures,  and  service  conditions.  To  test  the  adequacy  of  K{  and  flow  stress  criterion,  test 
data  from  UKAEA  (Tables  4  &  5)  and  GE  (Table  6)  have  been  used  to  compute  Kc,  <t\  the  results 
are  shown  in  Figs.  22  and  25.  If  the  pseudo-toughness,  Kc ,  can  be  considered  as  a  criterion,  it  should 
be  a  constant  for  different  flaw  lengths  (geometry  independent).  Similarly,  if  the  flow  stress  criterion  is 
workable,  a  close  correlation  should  exist  between  the  flow  to  failure  stress  ratio,  <r  */«•,,,  and  the  M 
factor  (Eq.  4).  Unfortunately,  the  expected  behavior  is  not  apparent  in  these  figures. 


For  the  case  of  surface  flaws,  complexities  similar  to  those  for  through-wall  flaws  also  exist.  How¬ 
ever,  it  has  been  pointed  out  by  Kiefner  [5]  that  surface  flaws  have  much  lower  transition  temperature 
than  through-wall  flaws.  This  is  a  reasonable  observation  because  the  ligament  is  subjected  to  much 
less  constraint  than  a  full  thickness  section.  Consequently,  for  the  same  material  operating  at  the  same 
temperature,  surface  flaws  are  more  likely  to  undergo  ductile  rupture  rather  than  a  through-wall  flaw. 
Unfortunately  evidence  on  testing  of  A 106  B  pipes  [34]  having  inner  or  outer  surface  flaws  does  not 
substantiate  this  suggestion  (Tables  7,  8).  The  flow  stress,  computed  from  the  equation: 


yield 


+  <T, 


2.4 


(15) 


ranges  from  44  ksi  to  45  ksi.  If  the  surface  flaws  behave  in  ductile  manner,  Mp  a- h  values  on  these  two 
tables  should  be  within  this  range  in  order  for  the  flow  stress  theory  (Eq.  4)  to  hold.  Figures  26  and  27 
summarize  this  comparison  and  two  observations  are  noted.  First,  the  scattered  data  points  indicate 
lack  of  correlation;  and  second,  the  definition  for  the  "flow  stress"  is  not  optimal.  If  the  constant  2.4  in 
the  above  equation  is  reduced,  i.e.  higher  flow  stress  is  computed,  better  correlation  can  be  attained.  In 
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0  4  8 

Flaw  Length  (inches) 

Fig  24— GE  axial  through-wall  flaw 
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comparing  through-wall  flaws  (Fig.  25)  with  surface  flaws  (Figs.  26  &  27),  it  is  interesting  to  note  that 
if  higher  flow  stress  value  are  used,  the  surface  flaw  results  display  better  correlation.  This  may  indicate 
the  surface  flaw  have  higher  ductility,  lower  transition  temperature,  as  discussed  previously. 

3. 3. 2. 2.  Circumferential  Cracks 

When  a  material  is  operating  on  the  shelf,  ductile  rupture  is  expected.  The  Battelle  [38]  test  on 
24  x  0.75  in.  A106  B  pipes  (Table  9)  suggest  the  same  limit  load  analysis  as  that  used  for  high  tough¬ 
ness  materials  is  applicable.  In  this  analysis,  the  load  carrying  capability  of  the  flawed  pipe  is  restricted 
to  the  limit  load.  The  applied  pressure  and  bending  loads  are  corrected  to  the  new  neutral  axis  position 
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due  to  the  existence  of  a  circumferential  flaw.  The  limit  load  can  be  computed  from  the  flow  stress  of 
the  material.  Figure  28  is  a  cross  plot  of  results  from  Table  9,  where  both  bending  moment  and  failure 
stress  are  presented  with  flaw  length.  It  is  noted  that  maximum  bending  moment  capability  does 
correspond  well  with  minimum  failure  pressure.  This  suggests  the  failure  description  of  circumferen¬ 
tial  flawed  pipes  is  quantitatively  correct. 
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Fig.  28— Failure  pressure  and  calculated  bending  moments  for  various  circumferential  surface  flaw  length 
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When  the  material  is  operating  at  room  temperature,  similar  analysis  technique  may  be  employed. 
Test  results  (45,  46]  on  6  in.  A106B  schedule  80  { ipes  are  tabulated  in  Tables  10  and  11.  To  further 
investigate  the  adequacy  of  the  limit  load  theory,  a  limit  moment  is  computed  for  0.0  and  1000.0  psig 
internal  pressure  at  various  through  crack  lengths  (2u)  Solid  lines  in  Fig.  29  represents  the  computed 
values  using  Eq.  11,  and  test  points,  represented  by  circles,  are  from  Tables  10  and  11.  A  fairly  good 
agreement  is  noted;  however  when  compared  with  the  high  toughness  material  (304  s.s)  predictions 
(Fig.  16),  the  experimental  result  displays  a  larger  deviation  from  theory. 

3.4  Recent  Advances  in  Elastic/Plastic  Fracture  Mechanics 

The  pseudo-toughness  criterion,  Kt,  was  developed  from  LEFM  theory,  modified  by  a  large  body 
of  experimental  data  and  it  has  been  applied  with  very  limited  success  to  pipe  rupture  studies  on 
medium  toughness  materials  operating  in  the  brittle  to  ductile  transition  region.  Previous  discussion 
has  also  indicated  that  the  flow  stress  theory  may  be  utilized  to  compute  fracture  behavior  of  piping 
made  of  high  toughness  materials  or  medium  toughness  materials  operating  at  shelf  temperature.  The 
desirable  situation  in  a  piping  integrity  investigation  is  to  have  material  criteria  that  are  a  counterpart  of 
Ku  for  brittle  fracture,  and  are  capable  of  characterizing  material  at  the  shelf  temperature  as  well  as  in 
transition  region. 

The  importance  of  recent  advancement  in  elastic-plastic  fracture  mechanic  is  not  only  to  provide  a 
better  understanding  of  the  fundamentals  of  fracture  propagation  processes,  but  also  to  clarify  the  limi¬ 
tations  and  applicability  of  simple  criteria,  e.g.,  T,  COA  or  "final  stretch"  that  have  been  derived  under 
assumptions  of  J initiation  criterion  [52,53]. 
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TABLE  10.  BEHAVIOR  OF  CIRCUMFERENTIALLY  FLAWED  6-INCH  A106B  SCHEDULE 
80  PIPES  UNDER  BENDING  AND  PRESSURE  STRESSES 
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2  «  (DEGREES) 

Fig  29— Comparison  of  limit  moment  predictions  with  experimental  results— A106B  piping 
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3.4.1  Stable  Crack  Growth  Criteria  Based  on  the  J-Controiied  Growth 

To  treat  elastic/plastic  fracture  initiation  and  propagation  problems,  crack  opening  displacement 
(COD),  J  controlled  growth  and  energy  method  are  candidates  currently  under  intensive  investigation. 
Although  they  have  been  proposed  at  different  stages  in  the  development  of  elastic/plastic  fracture 
theory,  there  is  strong  correlation  among  them  [47,48,49.50,51).  Because  of  the  underlying  assumption 
in  tnese  theories,  extensive  research  effort  is  being  carried  on  tr  establish  the  conditions  under  which 
they  can  be  reliably  utilized.  To  date,  there  is  no  single  theory  mat  can  pass  close  scrutiny.  However, 
provided  certain  requirements  are  satisfied,  it  is  evident  from  the  vast  amount  of  available  information 
ihat  the  /-controlled  growth  method  can  be  employed  as  a  useful  tool  to  define  structural  stability 
against  elastic/plastic  fracture. 

Following  the  introduction  of  the  /-integral  [7],  Ju.  has  been  accepted  as  a  elastic/plastic  fracture 
initiation  criterion  [10,11,12,13]  and  the  existence  of  a  HRR  (Hutchinson-Rice-Rosengren)  field  [8,9] 
has  been  identified  as  the  necessary  condition  of  a  /-dominated  stress  field.  In  the  /-integral  derivation, 
it  is  assumed  that  no  unloading  occurs  (deformation  theory  instead  of  incremental  theory)  and 
infinitcsmal  deformation  (as  opposed  to  finite  strain)  governs  the  kinematics  relationship.  Within  the 
requirements  of  crack  initiation,  apparently  both  assumptions  are  acceptable.  However,  once  the  crack 
starts  to  propagate,  further  large  deformation  at  the  crack  tip  is  expected  and  most  of  all,  unloading  is 
also  experienced.  Research  efforts  have  been  directed  to  extend  /  to  a  governing  propagation  parameter 
in  spite  of  its  underlying  assumptions.  Some  efforts  have  been  proven  successful. 

During  crack  growth  there  is  some  elastic  unloading  in  the  wake  of  advancing  crack  tip  where  the 
strains  are  strongly  concentrated  and  where  distinctly  non-proportional  plastic  deformation  occurs  near 
the  crack  tip.  Since  the  /-integral  is  theoretically  based  on  the  deformation  theory  of  plasticity  it  does 
not  adequately  model  any  of  these  aspects  of  plastic  behavior.  Although,  strictly  speaking,  /  is  res¬ 
tricted  to  the  analysis  of  stationary  cracks,  a  rationale  is  given  by  Hutchinson  and  Paris  [53]  for  use  of  / 
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to  analyze  crack  growth  and  stability  under  conditions  which  are  called  ./-controlled  growth.  Assuming 
fully  plastic  situations  for  generally  strain  hardening  materials,  they  considered  the  HRR  (Hutchinson- 
Rice-Rosengren)  type  strain  field  based  on  the  deformation  theory  for  a  growing  crack,  which  is: 

n  n 

KnJ"+>  r  "+l€u(d)  (16) 

where, 


K, ,  =  a  constant 

n  —  strain  hardening  coefficient 


r,9  =  local  polar  coordinates  at  the  current  crack  tip. 


The  strain  increments  under  a  simultaneous  increase  in  /  and  crack  length  are  written  as 


dt„  =  K„J 


n  + 1 


rt  +  1 


n 


n  +  1 


dJ  . 
J  * 


—  <u  + 


where 


(17) 


”  r  cos  9  i,j  +  sin  9  ~  (18) 

n  +  1  ov 

da  =  increment  of  crack  length,  and  A  a  <<  R. 

Equation  (17)  indicates  that  the  loading  increments  will  be  proportional  («/«„  ~  e(/)  if  the  first  term 
dominates,  which  may  be  assured  when 


If  there  exists  an  annular  region 


da 

r 


dJ 

<  <-y. 


where 


D  <  <  r  <  R, 


D 


dJ_  J_ 
da  J 


(19) 


(20) 


(21) 
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and  R  is  a  characteristic  dimension  for  HRR  field,  then  plastic  loading  is  proportional  and  HRR  singu¬ 
larity  is  dominant  in  this  region.  Therefore,  deformation  theory  still  remains  as  valid  procedure  for  r  > 
>  D  and,  the  significance,  for  purpose  of  the  present  discussion,  is  that  J  uniquely  measures  or  physi¬ 
cally  controls  the  fields  specified  by  Eq.  20  for  a  growing  crack.  It  may  also  be  assumed  that  if  predom¬ 
inantly  proportional  loading  occurs  throughout  most  of  the  singularity  region,  it  will  occur  outside  ( r  > 
R)  this  region  as  well  [53],  Thus  the  requirements  for  the  ./-controlled  growth  are  established. 


If  plane  strain  conditions  are  assumed  to  be  present,  which  may  be  defined  by  the  size  require¬ 
ment  criterion. 


Size^  25  J/o-0,  (22) 

the  /-Resistance  curve  is  size  independent  and  may  be  reasonably  configuration  independent.  Paris,  et 

al.  [54,59]  proposed  a  Tearing  Modulus  to  characterize  a  materials'  stable  tearing  property. 


_  _E_  dl 

o’  q  da 


(23) 


where  E  is  Young’s  modulus  and  a  is  the  crack  length.  The  behaviors  of  a  structure  during  crack 
growth  Uapp  >  / / )  is  determined  by  the  equilibrium  condition 


■/app  *^mat 

and  the  stability  criterion, 

/'app  <  T’ma.  =  stable  (24) 

/"app  >  Tmai  =  unstable 

where 


E  dJ  mai 

‘  mat  “  2  A*, 

ortf  da 

The  onset  of  instability  holds  when  rapp  -  Tmal. 


(25) 
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For  a  bi-linear  R-curve,  dJ / da  is  simply  the  slope  of  the  curve  during  crack  extension;  however, 
for  arbitrary  (non-linear)  resistance  curves,  it  should  be  noted  that  the  Tearing  Modulus  ( T)  is  depen¬ 
dent  on  the  instantaneous  J  value. 

Paris  and  his  colleagues  have  demonstarated  the  applicability  of  the  Tearing  Instability  Criterion 
for  center  cracked  tension  panels  [56],  HSST  intermediate  test  vessels  [56),  and  BWR  pipe  rupture  test 
[55],  The  effects  of  material  strain  hardening,  small  scale  yield  [54],  and  the  correlation  between  tear¬ 
ing  instability  and  Turner’s  t)  factor  [57]  have  also  been  investigated.  An  experimental  program  utiliz¬ 
ing  three-point  bend  specimens  was  recently  conducted  and  a  good  correlation  between  observed 
behavior  and  T-predicted  instability  response  was  noted  [54], 


Based  on  Rice’s  [601  Jx  flow  theory  for  an  ideally  plastic  material,  similar  studies  have  been 
conducted  by  Shih  et.  al.  [61].  They  considered  the  rate  of  change  of  the  strain  field  crack  growth  in 
the  form: 


da  r  da  ’  E  r 


R{0) 


SijW 


where 


(26) 


8  =  crack  opening  displacement 

R  (0)  =  a  measure  of  the  plastic  zone  size 
A  a  =  crack  growth;  <  <  R 

a-Q  =  flow  stress 

E  =  Young’s  modulus 

/y(0),gy(0)=  dimensionless  functions  of  order  unity. 
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The  first  term  in  Eq.  26  represents  the  strain  increment  rate  associated  with  crack  tip  blunting  while  the 
second  term  represents  the  strain  increment  rate  caused  by  the  crack  advance.  In  Eq.  26,  the  strains  at 

the  crack  tip  may  be  characterized  by  the  crack  opening  angle,  if  the  first  term  dominates  that  is, 

da 


db  o'o  . 

~r  >  >  ~pr  tn 
da  E 


R  (6) 


(27) 


Therefore,  Eq.  27  becomes  the  requirement  for  a  crack  opening  angle-controlled  crack  growth.  It 
should  be  noted  that  the  strain  fields,  Eq,  21  and  Eq.  26,  derived  from  two  different  approaches,  viz,  J j 
deformation  theory  and  J2  flow  theory,  have  a  very  similar  structure.  Therefore,  if  the  HRR  field 


grows  more  rapidly  than  the  advancing  crack,  either  — —  or  may  characterize  the  crack  tip  environ- 

da  da 

ment  for  a  growing  crack.  Based  on  the  crack  opening  angle,  Shih,  et  al.  [61],  proposed  a  tearing 
modulus. 


_£ 

a  o  da 


(28) 


as  a  parameter  to  characterize  the  stable  crack  growth  and  its  stability.  Extensive  studies  on  verification 
of  the  validity  of  the  crack  opening  angle  to  characterize  crack  growth  have  been  done  both  experimen¬ 
tally  and  numerically  by  Shih,  et  al  [61).  and  Hahn,  et  al.  [62] . 


Sorensen  [50]  and  Rice  and  Sorensen  [51]  formulated  a  ductile  crack  growth  criterion  which 
describes  the  critical  magnitude  of  a  crack  tip  opening  angle  during  the  crack  growth.  They  assumed  a 
Prandtl  slip-line  stress  field  centered  at  a  moving  crack-tip  in  an  elastic  perfectly-plastic  von  Mises 
material.  As  a  consequence  of  having  the  Prandtl  field  zone  translate  through  the  material  with  the 
advancing  crack,  they  derived  a  crack  tip  opening  angle  during  crack  growth  as. 


ft  ^  a  dJ  <r„  eR 

Au  <T„  da  P  E  A  a 


(29) 
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where 


cr„,£  =  as  defined  previously 

R  =  plastic  zone  size;  small  scale  yielding  is  assumed 
A  a  —  increments  of  crack  extension;  Aa  <  <  R  is  assumed 

8  =  crack  surface  displacement  at  the  previous  crack  tip 


a,fi  =  constants. 

If  a  stationary  crack  under  monotonic  loading  is  considered  (i.e.,  A  a  =  0),  Eq.  (29)  integrates  to  the 
crack  tip  opening  displacement, 

8 /  (30) 

<To 

The  material  dependent  non-dimensional  Constanta  correlates  8,  with  J  at  the  loading  stage  prior  to  the 
crack  growth  initiation.  Another  material  constant  /3  is  theoretically  defined  as  f3  =  4  (2—  v)/V3, 
where  v  is  Poisson’s  ratio.  However,  the  theoretical  value  of  /3  may  not  accurately  fit  the  relation  for 
the  results  obtained  numerically,  because  the  discretized  numerical  procedure  can  not  precisely  simulate 
the  continuous  crack  extension.  As  demonstrated  by  Sorensen  [50],  /3  can  be  reestimated  from  finite 

element  solutions  according  to  the  following  analysis. 


Let  the  crack  be  incremented  by  A  a  at  a  constant  external  load  level  so  that  the  first  term  in  the 
right  hand  side  of  Eq.  (29)  is  zero.  Assuming  the  plastic  zone  size,  R,  is  proportional  to  EJ/cr „  as  is  in 
the  case  of  monotonic  loading  of  a  stationary  crack,  (3  may  be  estimated  by 


8  E  a  , 

- - =  13  In 

A  a  (r  „ 


XcEJ 


(Tn&a 


where  X  is  a  proportionality  used  for  R,  and  may  be  fitted  together  with/3. 


(31) 
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Taking  S/A  a  to  be  a  material  dependent  constant  during  the  crack  growth,  and  assuming  small 
scale  yielding  (R  =  kEJ/al)  Sorensen  derived  the  following  crack  growth  criterion  from  Eq.  (29): 


dJ  f3  <*  o 


da 


In 


(32) 


dJ 


where  /«  denotes  J  at  steady  state  conditions  —  =  0,  so  that 

da 


.  <rfaa 
J„  =  — — —  exp 


keE 


E 

8 

P<T  o 

A  a 

(33) 


The  criterion,  Eq.  (32),  describes  the  J  requirement  necessary  to  continue  stable  crack  growth.  Equa¬ 
tion  (32)  is  integrated  to  obtain  the  current  crack  length  as  a  function  of  the  required  J  level,  i.e., 


(To  a  ~  a0 


4 

Ei  In 

ss 

~  Ei  In 

J 

he 

(34) 


where 


-oo  e~u 

Ei(x)  =  — —  tfw  (Exponential  integral) 

The  crack  growth  criterion  may  be  restated  such  that  the  necessary  and  sufficient  condition  to  sustain 
crack  extension  is  that  the  applied  J  equals  the  J  required  to  meet  the  growth  criterion  (Eq.  (34)),  thus 

JA(Q,a)  —  J  (a  —  a„)  (35) 

where  Q  denotes  applied  load.  Likewise,  instability  will  occur  when  applied  value  of  dJ/da  equals  or 
exceeds  the  dJ/da  required  to  meet  the  fracture  criterion,  thus 

dJA(Q,a)  ^  dJ(a  —  a0)  ,,,, 

7  ^ 

da  da 

Wnuk  (63]  also  formulated  an  instability  criterion  for  a  growing  crack  based  on  the  concept  of 
final  stretch,  which  is  defined  as  an  incremental  displacement  continually  generated  at  the  crack  tip  dur¬ 
ing  the  stable  crack  growth,  linked  with  the  Dugdale-Bilby-Cottrell-Swinden  model.  He  obtained  a 
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resistance  curve  for  plastic  zone  size  R  and  J  for  a  cracked  panel  of  infinite  width  that  experienced  a 
large  scale  yielding: 


or 


R  =  Rj 


K  ,+  l 


exp 


-i'"2 


j-j, 


a 

8  <Ty 

+  - 7Ta 

Xfln 

a 

-  -j  In2 

a 

<*o 

itE 

a0 

4 

a°\ 

(37) 


where  cry  is  yield  strength  and  the  subscript  /  implies  the  value  of  the  parameters  at  the  moment  of 
crack  growth  initiation.  The  parameter  A  measures  the  distance  of  any  state  during  the  stable  crack 
growth  from  the  ultimate  loss  of  stability.  Thus,  it  is  called  "stability  index”.  At  crack  growth  initia¬ 
tion,  the  index  A  is  defined  as, 


—  n E 

1  8cr  y 


dJ 


da 


9  Ja 


da 


(38) 


dJ  9/4 

where  —  is  a  measure  of  material  resistance  and  the  — —  is  an  applied  value.  It  is  obvious  that  A,  can 
da  da 

be  expressed  in  terms  of  tearing  modulus  T.  Stable  crack  growth  will  occur  if  A,  is  greater  than  zero, 
and  the  crack  will  become  unstable  at  the  point  where  the  A  drops  to  zero.  However,  for  some  types  of 
cracked  ductile  specimens  the  crack  may  not  grow  to  the  critical  size  of  instability.  In  such  a  case  A  will 
never  drop  to  zero  and  the  specimen  fails  through  a  ductile  tearing  mechanism.  Wnuk  combined  Eq. 
(37)  with  Eqs.  (35)  and  (36)  and  derived  the  following  expression  for  crack  length  at  the  terminal  ins¬ 
tability  for  a  central  crack  panel: 


where 


af=  a0exp[2A,] 
af  =  A  a  exp  [2(C-  1)1 


(39) 
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where  8  denotes  the  final  stretch. 


C=T~T^-\l"  (2e) 

A  a  8a- f  2 

The  terminal  toughness  is  given  by 


J,  =  J/K' 


1  +  flLMX2 
J,  nE  K‘ 


or 


j  - 

*  rt 


a„ 

1  +  ~- 


8o-  v 


Jj  IT  E 


(C- l)2-  (C-  \)ln 


a0 

+  -  In2 

a0 

A  a 

4 

A  a 

(40) 


The  advantage  in  these  equations  is  that  they  can  be  used  to  predict  the  critical  conditions  in  terms  of 
known  constants.  It  is  interesting  to  note  that  similarities  exist  in  the  mathematical  forms  between  the 
results  obtained  independently  by  Rice  and  Sorensen  and  by  Wnuk.  If  A  a  is  assumed  to  be  small  com¬ 
pared  to  the  plastic  zone  size  R ,  Wnuk’s  equations  can  be  modified  and  matched  to  those  by  Rice  and 
Sorensen. 


To  summarize,  although  based  on  two  different  theories  of  plasticity,  i.eN,  deformation  theory  and 
flow  theory,  the  requirements  of  ./-controlled  crack  growth  and  crack-tip  opening  displacement- 
controlled  growth  have  been  obtained  by  Hutchinson  and  Paris  [54]  and  Shih  et.  al.,  [611  respectively. 
If  these  requirements  are  met,  J  and/or  8  remain  as  valid  tools  to  characterize  stably  growing  cracks. 
Then  Eqs.  (24),  (36),  and  (38),  which  are  closely  related  to  one  another,  give  a  fracture  criterion  in 

terms  of  tearing  modulus  F,  /-resistance  — ,  and  stability  index  \  respectively.  In  addition,  the  insta- 

aa 

bility  condition  is  expressed  in  terms  of  terminal  crack  length  by  Eq.  (39),  and  in  terms  of  critical  J 
value  by  Eq.  (40). 

3.4.2  Stable  Crack  Growth  Criteria  Based  on  the  Griffith  Energy  Balance  Equation 

There  exist  other  types  of  studies  on  stable  crack  growth  whose  associated  failure  criteria  are  not 
directly  related  to  the  /-controlled  growth.  Some  of  these  are  based  on  the  Griffith  energy  balance  con¬ 
cept  extended  to  elastic-plastic  materials  undergoing  slow  crack  growth  in  finite  step  sizes,  A  a. 
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Neglecting  dynamic  and  thermal  effects  on  the  structure,  the  energy  balance  for  the  structure  dur¬ 
ing  a  step  of  crack  advance  may  be  written  as 

A  Wf  A  Wt  +  A  W.  +  A  W, 

~~  =  - - - 7  - - (41) 

A  a  A  a 

where 

A  Wf  =  the  external  work  increment  applied  to  the  structure 
A  We  =  the  elastic  energy  increment  stored  in  the  structure 
A  Wp  =  the  plastically  dissipated  energy  increment  in  the  structure 

and 


A  Wc  =  the  work  dissipated  due  to  separation  of  the  crack  surface  over  A  a. 


Postulating  A  a  is  finite,  define 


G*  = 


A  W}  -  A  Wc  -  A  Wp 
A  a 


and 


(42) 


= 


A  Wc 
A  a 


(43) 


then  G*  may  be  interpreted  as  the  rate  of  energy  available  for  creating  new  crack  surfaces,  and  G A  is 
the  rate  of  work  required  to  quasi-statistically  release  the  cohesive  tractions  holding  the  crack  surfaces. 
During  the  process  of  stable  crack  growth,  <7*  and  GA  are  equivalent.  Atluri,  et.  al.,  [641  and  Hahn,  et. 
al.,  [62]  studied  the  G*  by  finite  element  analysis  which  simulated  a  resistance  curve  from  existing 
experimental  data.  On  the  other  hand,  the  G *  has  been  studied  by  Hahn  et.  al.,  [62]  and  Kfouri  and 
Miller  [65].  As  predicted  by  Rice  [66]  that  at  least  for  non-hardening  materials,  both  G*  and 
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Another  candidate  for  a  criterion  to  characterize  stable  crack  growth  and  possibly  predict  the  onset 
of  fracture  is  a  crack  tip  nodal  force.  This  has  been  proposed  and  studied  by  Hahn,  et.  al.,  [62], 
Although  this  is  intuitively  appealing,  the  theoretical  basis  for  use  of  the  crack  tip  nodal  force  has  not 
yet  been  established. 

4.0  LEAK  BEFORE  BREAK  CRITERIA 

It  is  very  desirable  that  pipe  rupture  is  preceded  by  leakage  so  that  necessary  measures  can  be 
undertaken  to  remedy  the  situation.  Therefore,  there  are  two  aspects  to  be  considered:  (1)  the  possibil¬ 
ity  of  leak  before  break,  and  (2)  the  acceptable  leakage  rate.  The  state  of  art  of  pipe  rupture  in  the 
elastic/plastic  region  is  only  at  the  semi-empirical  stage.  Therefore,  effort  has  been  devoted  to  the 
definition  of  critical  crack  length  with  only  little  emphasis  to  the  mouth  opening,  which  is  one  of  the 
means  to  measure  leakage  rate. 

For  pipes  fabricated  from  austenitic  stainless  steel  or  ferritic  steels  (at  shelf  temperature),  pipe 
fracture  is  unlikely  to  occur  unless  the  structure  has  reached  its  limit  load  capability.  In  that  case,  duc¬ 
tile  tearing  is  the  failure  mode.  From  all  the  experimental  and  analytical  results  reviewed,  the  crack 
length  that  corresponds  to  limit  load  (flow  stress  criterion)  is  generally  very  large  and  excessive  leakage 
would  have  been  developed  well  before  crack  of  such  length  is  developed  [3,  29], 

For  pipes  fabricated  from  ferritic  steels  that  operate  in  the  transition  temperature  region,  mixed 
mode  failure  consisting  of  brittle  fracture  and  ductile  tearing  is  observed.  The  leakage  rate  is  deter¬ 
mined  by  the  degree  of  mouth  opening.  Analytically,  the  crack  opening  can  be  approximated  as  the 
sums  of  elastic  displacement,  plastic  correction  and  bulging  to  pipe  geometry.  The  following  derivation 
is  from  Ref.  29. 
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A  =  Mouth  opening  area  =  A\  (Elastic)  +  A2(PIastic)  +  A3  (Bulging),  (47) 

,  ,,  2(1  -  v2)  l~c  „ 

where  A ,  =  - - -  •■»/  —  Knc 

E  V  it 

4<r..c2 

A  2  =  — = —  In  Sec 
t 

and 

a}  = 

in  which, 

2c  =  crack  length 

v  «=  Poisson’s  Ratio 

£  =  Modulus 

c rv  =  yield  strength 

tr  —  applied  stress 

=  amount  of  bulging  (Ref.  29) 

K  —  stress  intensity  factor 

R  =  pipe  radius 

5.0  CONCLUSIONS 

An  extensive  literature  review  concerning  experimental  data  and  predictive  capability  of 
elastic/plastic  pipe  rupture  has  been  completed.  The  materials  used  for  nuclear  reactor  piping  can  be 
classified  as  high  toughness  and  medium  to  high  toughness  materials.  High  toughness  material,  includ¬ 
ing  austenitic  stainless  steel,  exhibits  no  brittle  to  ductile  transition  temperature  and  ductile  tearing  is 
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the  failure  mode.  Medium  to  high  toughness  materials,  including  ferritic  steels,  SA  106,  SA  333  and 
SA  561,  exhibit  high  ductility  at  the  shelf  temperature  and  fail  in  a  mixed  (brittle/ductile)  fracture 
mode  in  the  transition  temperature  region. 

The  predictive  tools  reviewed  are  classified  into  three  major  categories;  the  overall  structural 
response  computation,  the  semi-empirical  method,  and  rigorous  elastic/ plastic  fracture  mechanics. 
Lack  of  adequate  elastic/plastic  fracture  criteria  and  excessive  computational  expense  limit  the  applica¬ 
bility  of  overall  structural  response  computation.  For  the  cases  studied,  the  flow  stress  theory,  which 
assumes  pipe  rupture  from  limit  load,  is  adequate  to  predict  failure  of  high  toughness  materials  and 
medium  to  high  toughness  materials  at  the  shelf  temperature.  The  pseudo-toughness  approach,  Kc,  has 
been  used  to  describe  pipe  failure  in  the  transition  temperature  region  with  limited  success.  Therefore, 
for  those  materials  which  operate  in  the  transition  region,  the  real  solution  is  to  ultilize  rigorous 
elastic/plastic  fracture  mechanics  methodology.  Among  the  current  theoretical  developments  in  the 
analysis  of  advancing  cracks,  the  ./-controlled  growth  methods  are  one  of  the  most  promising  and  most 
widely  accepted.  All  the  criteria  proposed  to  data,  i.e..  Tearing  Modulus,  COA,  final  stretch,  and 
Stability  Index,  are  capable  of  describing  a  small  amount  of  crack  growth  and  they  are  all  mathemati¬ 
cally  related.  To  permit  quantitative  assessment  of  safety  factors,  ./-integral  resistance  curves  of  current 
nuclear  piping  materials  should  be  generated. 

The  high  fracture  toughness  levels  associated  with  existing  nuclear  piping  configurations  and 
materials  suggests  that  extensive  plastic  deformation  will  precede  failure.  No  sudden  (catastrophic) 
ruptures  are  expected  and  a  favorable  leak-before-break  condition  should  prevail.  The  mouth  opening 
on  the  pipes  of  these  materials  in  the  transition  region  can  be  approximated  by  linear  elastic  mechanics 
modified  by  plasticity  and  geometry  corrections. 
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TABLE 

Summay  of  flat-plate  crack-extension  measurements  performed  by  Kthara  et  at.  on 
19.7-in  -long  by  15.1-tn  -wide  hot-rolled  steel  plates  with  X 14 -in.-long  edge  cracks. 


Plate  thickness,  in. 

Test  temp,  C 

e*.  psi 

ksi 

0.25 

-118 

20,300 

73 

-138 

18.000 

65 

-158 

14,500 

52 

-170 

15,600 

56 

-196 

8,170 

29 

0.375 

-118 

25,200 

90 

-138 

20,200 

73 

-158 

15,600 

56 

-180 

12,400 

45 

0500 

-118 

28,400 

100 

-138 

23,800 

85 

-158 

17,000 

61 

-180 

9,900 

36 

1,1  Kt=  1.14  er*(jrc)*,  no  plasticity  correction  was  necessary  since  o'! a,  <  0.4  in  all  cases. 


Compilations  of  pressure-vessel  test  data 


table 

Data  summary  for  Anderson  and  Su/fica/i  for  aluminum-olio y  vessels 


Mater  al  properties 

Vessel  geometry 

Test  conditions 

Calculated  values' 

M 

Designation 

or,  ksi 

ov,  ksi 

R.  in. 

t,  in. 

Rft 

Temp,  C 

c.  in. 

*S. 

c’/Rl 

9>> 

(*;*»«:*,)- 1  X 

10- *•(»».»)  (IbT*) 

Aluminum  20I4-T6 

68 

79 

281 

0.06 

46.8 

Room  Temp. 

006 

644 

0.02 

142/1  55 

5.29/  8  26 

Ditto 

68 

79 

2.81 

006 

46.8 

Ditto 

025 

340 

0.37 

1.2  / 1.14 

9  17/  9  66 

Ditto 

68 

79 

2.81 

006 

46.8 

Ditto 

050 

20.6 

1.48 

1.13/1.09 

13.40/13.76 

Ditto 

68 

79 

181 

006 

46.8 

Ditto 

1.00 

9.7 

5.93 

1.08/1.06 

31  32/3191 

Ditto 

68 

79 

2  81 

0.06 

46.8 

-196 

005 

71.6 

0014 

1  75/142 

7.09/  8.74 

Ditto 

82 

93.9 

2.81 

0.06 

46  8 

-196 

0.07 

70.6 

0029 

1.74/1.41 

5.24/  6.47 

Ditto 

82 

93.9 

2.81 

0.06 

46.8 

-196 

010 

63.7 

0059 

1.49/1. 3 1 

5.26/  5.98 

Ditto 

82 

93.9 

2.81 

006 

46.8 

-196 

012 

58.5 

0.085 

1.39/126 

5.57/  6.15 

Ditto 

82 

93.9 

2.81 

0.06 

46.8 

-196 

015 

512 

013 

1.29/1.20 

6.03/  6  49 

Ditto 

82 

93.9 

281 

0.06 

46.8 

-196 

0.20 

47.4 

023 

1 .25/1.17 

5  66/  6.03 

Ditto 

82 

93.9 

281 

006 

46.8 

-196 

0.25 

40.1 

037 

1.19/M  3 

665/  700 

Ditto 

82 

939 

2.81 

0,06 

46.8 

-196 

037 

30.2 

081 

I.14/1.IO 

8.27/  8.57 

Ditto 

82 

939 

2.81 

0.06 

46.8 

-196 

050 

231 

1.48 

Ml/108 

10.75/11.04 

Ditto 

82 

93.9 

2.81 

0.06 

46  8 

-196 

062 

186 

128 

1.09/107 

13.61/(3.86 

Ditto 

82 

93.9 

2.81 

0.06 

468 

-196 

087 

14.4 

4.50 

1.09/1.07 

16.10/16  48 

Ditto 

82 

93.9 

181 

0.06 

46.8 

-196 

1.00 

113 

5.80 

1.07/1.03 

23.30/23.68 

Ditto 

90  5 

939 

181 

0.06 

46.8 

-253 

005 

812 

0.014 

1.97/1.77 

4.78/  5.32 

Ditto 

90  5 

93.9 

181 

0.06 

46.8 

-253 

012 

63.4 

008 

1.38/1.34 

4.78/  4.92 

D..to 

90  5 

93.9 

181 

0.06 

46  8 

-253 

025 

396 

037 

1.16/ 1.14 

6.99/  7.12 

Ditto 

90  5 

93.9 

181 

006 

46  8 

-253 

0.37 

320 

0.81 

1.15/M4 

7.30/  7.37 

Ditto 

90  5 

93.9 

281 

006 

468 

-253 

050 

210 

1.48 

1.08/1  08 

13.36/13.38 

Ditto 

905 

93.9 

2.81 

006 

468 

-253 

0.63 

198 

135 

1.I1/M0 

II  60/11.70 

Ditto 

90  5 

93  9 

181 

006 

468 

-253 

087 

13.2 

4.52 

1.08/1.07 

1945/196) 

Ditto 

90  5 

939 

181 

006 

468 

-253 

100 

119 

593 

1.08/1  08 

20  80/20.80 

Two  values  are  computed  for  and  rrft1  rrr<p,  corresponding  to  the  upper  and  lower  bounds  for  d  :d  «er. 
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TABLE 

Data  summary  for  Getz,  Pitrce,  and  Cahet  jjr  aluminum-alloy  tessels 


Material  propel tics  Vessel  geometry  Test  conditions  Calculated  val vies'" 


Designation 

Of,  ksi 

Of,  ksi 

R,  in. 

i,  in. 

RJt 

Temp,  C 

c,  in. 

cJJ,  ksi 

<■*/*' 

Vs 

10  ,0(in.*)  (lb. '  * 

Aluminum  2014-T6  90.5 

93.9 

3.0 

006 

50 

-253 

0.05 

81 

0.013 

1.861.70 

5.21/  5.71 

Ditto 

905 

93.9 

3.0 

0.06 

50 

-253 

0.05 

83 

0.013 

1.97/1.80 

4.70/  5.14 

Ditto 

905 

93.9 

30 

006 

50 

-253 

0.05 

86 

0.013 

254/2.05 

3  39/  4  20 

Ditto 

9.05 

93.9 

3.0 

0.06 

50 

-253 

0.07 

73 

0.027 

1.53/1.46 

5.57/  5  84 

Ditto 

90  5 

93.9 

3.0 

006 

50 

-253 

0.12 

63 

0.08 

1.36/1.31 

491/  5  10 

Ditto 

90.5 

939 

3.0 

006 

50 

-253 

0.25 

40 

0.34 

1-14/1.13 

6.98/  7.04 

Ditto 

90.5 

93.9 

3.0 

0.06 

50 

-253 

0.25 

39 

0.34 

1.13/1.12 

7.41/  7.47 

Ditto 

90.5 

93.9 

3.0 

0.06 

50 

-253 

037 

33 

076 

1.12/1.1  J 

7.05/  7  18 

Ditto 

90.5 

939 

3.0 

0.06 

50 

-253 

0.50 

21 

1  38 

1.06/1  06 

13.63/13.63 

Ditto 

905 

93.9 

30 

006 

50 

-253 

0.61 

20 

2.06 

1  08/1.07 

12.08/12.20 

Ditto 

90.5 

93.9 

30 

0.06 

50 

-253 

0.87 

14 

4.20 

106/1.06 

17.63,17  63 

Ditto 

90.5 

93.9 

30 

0.06 

50 

-253 

1.0 

11 

5.50 

1.05/1.04 

25.06/25.18 

Ditto 

90.5 

939 

3.0 

0.06 

50 

-253 

1.0 

13 

5.50 

1.07/1.06 

17.63/17  79 

Ditto 

82.0 

93.9 

3.0 

0.06 

50 

-196 

0.06 

70 

0.02 

1.65/1.40 

6.5  /  7.7 

Ditto 

82.0 

93.9 

3.0 

0.06 

50 

-196 

006 

72 

0.02 

1.80/1.45 

5.69/  7.06 

Ditto 

820 

93.9 

3.0 

0.06 

50 

-196 

006 

75 

0.02 

210/1.50 

4.49/  6.29 

Ditto 

820 

93.9 

3.0 

0.06 

50 

-196 

0.07 

72 

0027 

1  90/1.45 

4.62/  6.05 

Ditto 

82.0 

93.9 

3.0 

0.06 

50 

-196 

0.10 

6$ 

0055 

1.50/1.35 

5.02/  5.58 

Ditto 

82.0 

939 

3.0 

0.06 

50 

-196 

0.14 

60 

0.10 

1.40/1.30 

4.51/  4  86 

Ditto 

820 

93.9 

3.0 

006 

50 

-196 

0.15 

54 

0.12 

1  35/1.20 

5  39/  6  06 

Ditto 

82.0 

93.9 

3.0 

0.06 

50 

-19f 

0.20 

45 

0.22 

1.20/1.15 

6.55/  6  u4 

Ditto 

820 

93.9 

3.0 

0.06 

50 

-196 

0.20 

50 

0  22 

1.30/1  20 

4.89,  5.30 

Ditto 

820 

93.9 

3.0 

0.06 

50 

-196 

0.25 

39 

0.34 

1.20  1.10 

6  98/  7.62 

Ditto 

82.0 

93.9 

3.0 

0.06 

50 

-196 

0.25 

41 

0.34 

1.20,1. 10 

6.31/  6  89 

Ditto 

82.0 

93.9 

30 

006 

50 

-196 

0.37 

30 

0.76 

1.10  l  10 

8.69-  8.69 

Ditto 

820 

93.9 

3.0 

0.06 

50 

-196 

0.50 

23 

1.38 

1.10/1.07 

10.95,11.26 

Ditto 

82.0 

9X9 

3.0 

0.06 

50 

-196 

100 

20 

5  50 

1.05  1.04 

30  33  50.60 

Ditto 

82.0 

9X9 

3.0 

0.06 

50 

-196 

1  00 

11 

5.50 

1.06  1  04 

24.89, 2537 

*“*  T»o  values  are  computed  for  <pj  and  ej*  iu^p,  corresponding  to  the  upper  and  lower  bounds  for  e  6  =  ot,  d  =  <xt 


TABLE 

Data  summary  for  Duffy  ei  ai  for  ductile  track  extension  in  steel  pipes 


Material  properties  Vessel  geometry  Test  conditions  Calculated  values'" 


Designation*** 

Of,  ksi 

Ov,  ksi 

A,  in. 

i,  in. 

R/t 

Temp,  C 

c.  in 

c1  Rt 

<Pj 

(ffJJ-Vtf,)'  *  X 

10  ‘°(itL,)(Ib'J 

*•- : 
lO'10 

Steel  RR.  TR,  BB 

60 

80 

15 

0.375 

40 

-20-24 

0.5 

70.6 

0044 

a, 2. 00 

0,0.638 

2.00 

Ditto 

60 

80 

15 

0.375 

40 

-20-24 

0.5 

698 

0  044 

i/200 

0/0.654 

205 

Ditto 

60 

80 

15 

0.375 

40 

-20-24 

165 

56.2 

0.483 

x/1.85 

0/0.331 

3  16 

Ditto 

60 

80 

15 

0.375 

40 

-20-24 

1.65 

55.8 

0483 

X/l  85 

0/0335 

320 

Ditto 

60 

80 

15 

0.375 

40 

-20-24 

225 

46.  B 

090 

x/1.70 

0/0.380 

4  50 

Ditto 

60 

80 

15 

0.375 

40 

-20-24 

2.70 

42.6 

1.30 

x/l  85 

0,t>.3Sl 

*  5.50 

Ditto 

60 

60 

IS 

0375 

40 

-20-24 

3.20 

388 

1  82 

x,20 

0,0329 

6.64 

Ditto 

60 

80 

15 

0.375 

40 

-20-24 

4.40 

282 

3.34 

x/l  60 

0/0.569 

12.50 

Ditto 

60 

80 

IS 

0.375 

40 

-20-24 

4.40 

278 

334 

oc/I.SS 

0/0604 

1290 

Dilto 

60 

80 

15 

0375 

40 

-20-24 

440 

27.6 

3.34 

x  / 1.55 

0/0613 

13.10 

Ditto 

60 

80 

15 

0.373 

40 

-20-24 

4.40 

27.6 

3.34 

x/l. 55 

0/0.613 

13  10 

Steel  AF 

68 

64 

15 

0.375 

40 

-16-16 

2.65 

48.7 

1.25 

x/2.20 

0/0230 

4.22 

Ditto 

66 

64 

15 

0.375 

40 

-16-16 

2.65 

48.1 

125 

x/2.10 

0/0247 

4.32 

Ditio 

66 

64 

15 

0.375 

40 

-16-16 

2.65 

474 

1.25 

x/1.90 

0/0.281 

4  45 

Ditto 

66 

84 

15 

0.375 

40 

-16-16 

440 

316 

3.34 

x/l. 75 

0/0414 

10.00 

Ditto 

68 

84 

15 

0.375 

40 

— 16-16 

7.50 

165 

100 

19  /l. 35 

0.82  /l  156 

37.20 

Ditto 

66 

64 

IS 

0.375 

40 

-16-16 

10.00 

1 1.4 

17  6 

I  59/1.25 

1.54  /|. 96 

7690 

Steel  AC,  AD 

53 

75 

15 

0.375 

40 

2-23 

35 

330 

2.18 

x/l  .65 

0/0.506 

9.18 

Ditto 

53 

75 

15 

0375 

40 

2-23 

3  5 

328 

2.18 

x/l  65 

0,0  512 

9.30 

Ditto 

53 

75 

15 

0.375 

40 

2-23 

3.5 

32.4 

218 

x/l  62 

0/0.535 

9.53 

Steel  UU 

61.7 

70 

306 

025 

1225 

-2-0 

2.20 

25.3 

635 

x/x 

0 

1S60 

Ditto 

61.7 

70 

306 

0.25 

1225 

-2-0 

3.55 

15  3 

164 

x/x 

0 

42.70 

Steel  GP 

51 

73 

13 

0.281 

46.4 

17 

255 

37.9 

1  78 

x/x 

0 

690 

Steel  AH 

60 

80 

16 

0.406 

45 

-8-4 

2.7 

439 

10 

x/l  55 

00  394 

5  10 

Ditto 

60 

80 

16 

0.406 

45 

-  8  -4 

2.7 

44.7 

10 

x/l  57 

00376 

500 

Dilto 

60 

80 

16 

0.406 

45 

-8-4 

2.7 

46.0 

1.0 

X/l  65 

0/0337 

470 

Steel  VY 

62 

74 

18 

0  861 

21 

62 

485 

336 

1  50 

26  ,1  5 

0  223  0  387 

880 

Two  values  are  computed  for  <p,  and  e£J  *c<p,  corresponding  to  the  upper  and  lower  bounds  for  do*  «T.  «r  “  <rt.. 
^  The  steels  employed  in  this  study  are  X-50  and  X-60  grade  line-pipe  steels. 
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TABLE 

Data  summary  for  Nichols  et  al  for  ductile  and  semibrittle  crack  extension  tn  steel,  pressure  vessels 

Material  properties  Vessel  geometry  Test  conditions  Calculated  values'*’ 


Designation 

a,,  ksi 

oL.  ksi 

R.  in. 

r,  in. 

R;t 

Temp,  C 

r,  in 

<rj.  ksi 

c\'Ri 

*3 

1  * 

10  ,0(m.,)(lb  •) 

#•  *  a 

10  ,#l«n*(lb. 

0  36  C  Steel1*1 

34  5 

69.5 

30 

L0 

30 

1-51 

3 

27.6 

0.3 

13  /I  12 

0.6  /I  24 

Ditto 

34  5 

69.5 

30 

1.0 

30 

1  51 

3 

312 

0.3 

oc/1.16 

0/0.88 

Ditto 

345 

69  5 

30 

10 

30 

1-51 

6 

188 

1.2 

1.8  /1.09 

0.83/1  37 

Ditto 

34.  J 

69.5 

30 

1.0 

30 

1  51 

6 

17.9 

1.2 

161 1  08 

1.02/1.53 

Ditto 

345 

69.5 

30 

10 

30 

1-51 

6 

21.0 

1.2 

at/1.12 

0/1.07 

Ditto 

345 

69  5 

30 

10 

30 

1-51 

6 

15  9 

1.2 

1.4  /1. 06 

1.49/1.98 

Ditto 

34.5 

69.5 

30 

1.0 

30 

1-51 

12.37 

96 

5.0 

l  5  /1.07 

1.86/161 

Ditto 

34.5 

69.5 

30 

1.0 

30 

62-88 

3 

330 

a3 

— 

— 

9.8 

Ditto 

34.5 

69  5 

30 

to 

30 

62-88 

6 

25.7 

1.2 

— 

— 

ISO 

Ditto 

34.5 

69  5 

30 

1.0 

30 

62-88 

6 

27.7 

1.2 

— 

— 

13.0 

Ditto 

34  5 

69  5 

30 

1.0 

30 

62-88 

12 

117 

4.8 

— 

— 

610 

Ditto 

34.5 

69.5 

30 

1.0 

30 

62-88 

12 

15.2 

4.8 

— 

— 

43  28 

Ditto 

34  5 

69.5 

IB 

10 

18 

10-50 

6 

13.9 

2.0 

1.4  /1.06 

1.96/2  59 

Ditto 

34.5 

69  5 

18 

10 

18 

10-50 

6 

17.4 

2.0 

115/108 

081/1.62 

D-.tto 

34.5 

69.5 

18 

10 

18 

79 

6 

23.4 

10 

— 

— 

18.2 

Ditto 

34.5 

69  5 

57 

10 

57 

17 

6 

210 

0.63 

1.55/1.07 

a7  /1 .02 

Ditto 

34.5 

69.5 

57 

1.0 

57 

80 

6 

27.1 

0.63 

— 

— 

13.5 

0.13  C  Steel1®* 

40.0 

63  5 

30 

10 

30 

16-79 

6 

29.9 

1.2 

— 

— 

11.1 

Ditto 

400 

63  5 

30 

10 

30 

16-79 

6 

31.8 

1.2 

— 

— 

9.88 

Ditto 

40.0 

635 

30 

1.0 

30 

16-79 

6 

29.4 

1.2 

— 

— 

1156 

Ditto 

40.0 

63.5 

30 

1.0 

30 

16-79 

6 

31.8 

1.2 

— 

— 

9.88 

Ditto 

400 

635 

'0 

10 

30 

16-79 

12 

19.5 

4.8 

— 

— 

262 

0  16  C  Steel’*' 

310 

64.0 

30 

10 

30 

39 

6 

28.5 

12 

— 

13.5 

T wo  values  are  compared  for  <p,  and  cr^:  ncipj  corresponding  to  the  upper  and  lower  bounds  for  a :  6  <rr.  6  =*  ou. 

P’-isn  carbon  steel;  C  0  36®;,  Mn  0  44.-0.4d'*..  St:  0  10-0  13®.  The  mode  of  crack  extension  in  this  steel  was  100  percent  ductile  shear  above  51  C.  semi* 

brittle  below  51  C 

Aluminum  grain-refined  steel;  C  0.13*..  Mn  -  1.14®..  Si  0  12®'..  Crack  extension  mode  was  100  percent  ductile  shear  in  all  cases. 

Stiicra-clled  steel;  C:  0.16®.,  Mn  1.22*,.  Si:  0.20%  Crack  extension  mode  was  100  percent  ductile  shear  in  all  cases 


TABLE 

Data  summary  for  Kihara.  Ikeda,  and  Iwanga  for  brittle-steel  vessels 


Material  properties  Vessel  geometry  Test  conditions  Calculated  values 


Designation 

or,  ksi 

<rr.  ksi 

R.  in. 

t,  in. 

R/t 

Temp,  C 

c,  in. 

*2.  ksi 

cV*i 

Steel'®* 

115 

125 

4.3 

0.25 

17 

— 196 

2.44 

3.98 

5.56 

1.00 

8264 

Ditto 

115 

125 

4.3 

025 

17 

-196 

196 

4.95 

3.59 

1.00 

66  22 

Ditto 

115 

125 

4.3 

0.25 

17 

-196 

US 

8.25 

1.70 

1.00 

346 

Ditto 

115 

125 

43 

0.25 

17 

-196 

086 

132 

0.69 

1.00 

21.17 

Ditto 

IIS 

125 

8.0 

0.25 

32 

-196 

3.22 

4.95 

5.18 

1.00 

4048 

Ditto 

115 

125 

8.0 

0.25 

32 

-196 

253 

5.60 

3.20 

1.00 

40.16 

Ditto 

115 

125 

80 

0.25 

32 

-196 

1.85 

8.52 

1.71 

1.00 

23  69 

Ditto 

115 

125 

80 

0.25 

32 

-196 

1.16 

960 

0.67 

1.00 

29.76 

Ditto 

115 

125 

6.0 

025 

25 

-196 

165 

7.56 

1.80 

1.00 

33.78 

Ditto 

IIS 

125 

64 

0375 

17 

-196 

2.0 

4.40 

1.66 

1.00 

81.96 

Ditto 

115 

125 

64 

0375 

17 

-196 

20 

5.70 

1.66 

1.00 

4901 

Ditto 

115 

125 

43 

0.375 

11.4 

-196 

165 

7  40 

1.68 

100 

35.21 

Ditto 

115 

125 

80 

0.375 

21.3 

-196 

224 

9.90 

1.67 

1.00 

1442 

Ditto 

115 

125 

4.3 

0.50 

8.6 

-196 

230 

3.42 

544 

1.00 

176.0 

Ditto 

115 

125 

43 

0  50 

86 

-196 

263 

3  50 

3.20 

1.00 

99.0 

Ditto 

115 

125 

4.3 

050 

86 

-196 

194 

4.20 

1.75 

1.00 

92.59 

Ditto 

115 

125 

4.3 

0  50 

86 

-  196 

1  16 

920 

0.62 

1.00 

3L46 

Ditto 

115 

125 

so 

050 

16 

-196 

3.62 

240 

3.27 

1.00 

153.8 

Ditto 

115 

125 

8.0 

0  50 

16 

-196 

2.33 

480 

1.60 

100 

54.64 

Ditto 

115 

125 

80 

0.50 

16 

-196 

1.55 

7.20 

0.60 

100 

39.58 

Ditto 

115 

125 

64 

0.50 

12.8 

-196 

234 

490 

1.71 

1.00 

56  49 

w  Hot-Rolled  Steel,  C: 

0.25%.  Si : 

0.02%.  Mn : 

085%. 
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TABLE 

Data  summary  for  AlmonJ  rt  at.  for  JuctUe-st**l  vessels 


Material  properties 

Of.  ksi 

Vessel  geometry 

Test  conditions 

°M-  ksi 

Calculated 

values 

iotfnctp,)- 1  x 

10  ‘“(m.’HIb.  ;) 

10  ,c 

Designation 

<rt.  ksi 

R.  in. 

«,  in. 

Rfl 

Temp,  C 

c,  in. 

c*fRi 

Steel4** 

45.7 

660 

2.5 

0.5 

50 

-25 

1.125 

32.0 

1 01  — 

_ 

9.75 

Ditto 

45.7 

660 

2.5 

05 

5.0 

-25 

1.125 

35.5 

1  01 

-• 

7.94 

Ditto 

45.7 

660 

2.5 

0.5 

5.0 

-25 

1.125 

360 

101 

-- 

7.70 

Ditto 

47.0 

70.6 

2.5 

0.5 

50 

-5  5 

1  125 

38.0 

1.01 

- 

694 

Ditto 

470 

706 

2.5 

05 

5.0 

-5-5 

1  125 

36  5 

101  - 

- 

7.50 

Ditto 

470 

70.6 

2.5 

0.5 

5.0 

-5-5 

1.125 

35.5 

1.01 

— 

7.94 

Ditto 

49 

77 

2.5 

0.5 

5.0 

-68 

1.125 

425 

1.01 

-- 

5.54 

Ditto 

66 

88 

2.5 

0.5 

50 

-120 

1.125 

45.7 

101 

— 

480 

M  Hot  Rolled  Steel;  C:  0.14%,  Si:  0.26%.  Mb: 0.47%.  Crack  extension  occurred  by  ductile  fibrous  mode,  at  least  initially. 


TABLE 

Daw  tmmrrujrv  for  Peters  and  Kuht  for  aimminum-itltoy  vessels 


Material  properties 


Vessel  geometry  Test  conditions 


Calculated  values 


Designation 

<rv,  ksi 

ksi 

K  in. 

t,  in. 

Rft 

Temp.  C 

c,  in. 

of,,  ksi 

50(t :  R:)tanh 
( R  500 

Vs 

(ol'-r.cvs)'  ’  * 

10'  l0(in.Jj(!b.  ‘ ;J 

Aluminum  2024-T3 

365 

63.0 

14.4 

0.015 

960 

Room  Temp 

0.31 

41  3 

0023 

x/l  25 

0  4.81 

Ditto 

36.5 

65.0 

14.4 

0015 

960 

Ditto 

064 

29  8 

0  (*98 

1  74  1)2 

3.21  4.99 

Ditto 

36.5 

650 

14.4 

0.015 

960 

Ditto 

1  25 

20  4 

ft  37 

l  ?:  !  07 

4  67,  5  7] 

Ditto 

365 

65.0 

14.4 

0.015 

960 

Ditto 

255 

It  3 

1  56 

118  1  W 

8.28  9.59 

Ditto 

36.5 

65.0 

14.4 

0.015 

960 

Ditto 

3.85 

82 

3  57 

1  18, 1  05 

104;  11/1 

Ditto 

36.5 

65.0 

36 

0025 

144 

Ditto 

050 

219 

0  96 

2.62, 1  1 5 

5  06  U  S 

Ditto 

365 

65.0 

36 

0.025 

144 

Ditto 

1.0 

10.4 

3  $4 

1.41  1.08 

20  85  27.23 

Ditto 

36.5 

650 

3.6 

0.025 

144 

Ditto 

2  81 

3.4 

29  50 

125  106 

7R.43  92  48 

Ditto 

36.5 

650 

3.6 

0.015 

230 

Ditto 

031 

29.3 

0  37 

x  1.18 

0  10  13 

Ditto 

36.5 

65.0 

3.6 

0.015 

230 

Ditto 

060 

172 

1  ?R 

1  6  ,1  11 

li  20  16 14 

Ditto 

365 

65.0 

3.6 

0015 

230 

Ditto 

1  20 

87 

5  50 

l  39,1  08 

25  24  52  4S 

Ditto 

36.5 

650 

3.6 

0.012 

300 

Ditto 

0.49 

229 

0  92 

x  1  16 

0  10  68 

Ditto 

36.5 

65.0 

3.6 

0.012 

300 

Ditto 

10 

128 

385 

2  18  i  14 

S  90  1  ■'.03 

Ditto 

36.5 

650 

36 

0.012 

300 

Ditto 

20 

5.7 

15.40 

1  46  1  09 

33  57  44  9 T 

Ditto 

365 

65.0 

3.6 

0012 

300 

Ditto 

1.90 

54 

13  90 

1.32  1.07 

43  5?  53  7 

Ditto 

36  5 

65.0 

36 

0012 

300 

Ditto 

048 

22.8 

0  88 

2.62  1  15 

4  8t  11  09 

Ditto 

36.5 

65.0 

3.6 

0.012 

300 

Ditto 

0.95 

10  2 

3*8 

1.34  107 

2399  30  yf 

Ditto 

365 

65.0 

36 

0012 

300 

Ditto 

028 

402 

0  05 

X/J  24 

0.  56" 

Ditto 

36.5 

65.0 

36 

0.012 

300 

Ditto 

012 

41.4 

0.05 

x,  1.26 

0  12.25 

Ditto 

365 

650 

3.6 

0012 

300 

Ditto 

j  25 

306 

0.25 

x  1.16 

0  11/; 

Ditto 

365 

65.0 

3.6 

0.012 

300 

Ditto 

0.25 

26.4 

0.25 

1.41/M  1 

12  96  16.4* 

Ditto 

36.5 

65.0 

3.6 

0.012 

300 

Ditto 

0.48 

20  0 

0  80 

1.57/1.10 

10  56  1  5  07 

Ditto 

36.5 

65.0 

3.6 

0.012 

300 

Ditto 

0.93 

10.7 

3  35 

1  37  1.08 

21/8  27  63 

Ditto 

36.5 

65.0 

36 

0.012 

300 

Ditto 

191 

4.9 

1250 

1.22,1  66 

56  92  65.51 

Ditto 

36.5 

65.0 

3.6 

0012 

300 

Ditto 

3.80 

26 

50.0 

1.24/1.66 

99  56  116  4* 

Ditto 

36.5 

63.0 

3.6 

0.012 

300 

Ditto 

0.48 

20.4 

080 

162/Ml 

9*2  14  3* 

Ditto 

36.5 

650 

36 

0012 

300 

Ditto 

0.48 

21.0 

0.80 

1.71/1.12 

8.^9/13.42 

Ditto 

36.5 

6S.0 

3.6 

006 

660 

Ditto 

0  12 

44.4 

0.05 

x/1.32 

0, 10  39 

Ditto 

36.5 

65.0 

36 

006 

600 

Ditto 

0.23 

27.3 

0.15 

1.56/1. 10 

11*9  16  86 

Ditto 

36.5 

65.0 

3.6 

0.06 

600 

Ditto 

048 

20.7 

0  80 

1.67/1.11 

9.26  13.94 

Ditto 

365 

65.0 

36 

0.06 

600 

Ditto 

0.98 

10.2 

3  50 

1  34/1.07 

23.32/29.2 

TABLE 

Data  mtmmary  for  Peters  mi  Kuhn  for  ahminmm~alloy  vtssrh 


Materials  properties  Vessel  geometry  Test  conditions  Calculated  values'*’ 


Designation 

«r.k« 

»i/.ksi 

A.  in. 

r,  in 

Rt 

Temp.  C 

c,  in. 

ksi 

50(r,/K,|  unh 

<Pi 

(»•’  «<*.)■'  * 

(A/501) 

10- "(in .‘)(H> 

Aluminum  7075-T6 

65 

80 

36 

0016 

225 

Room  Temp 

0.33 

19.4 

040 

106/104 

24  68/  24  65 

Ditto 

65 

to 

3.6 

0.016 

225 

Ditto 

065 

11.7 

160 

105/103 

3408,  34  74 

Ditto 

65 

80 

3.6 

0016 

225 

Ditto 

1  28 

5.5 

4  50 

1021  01 

80.60  81.43 

Ditto 

6$ 

80 

3.6 

0025 

144 

Ditto 

0  50 

166 

0  50 

105/103 

22  01  22  44 

Ditto 

65 

80 

36 

0.025 

144 

Ditto 

10 

8.4 

390 

106/1  03 

42.57/  43.81 

Ditto 

65 

80 

36 

0025 

144 

Ditto 

2.0 

3.7 

15  50 

1  04,1  02 

111  84  11403 

Ditto 

65 

80 

144 

0016 

960 

Ditto 

0.31 

372 

002 

2  02/1  39 

3.67/  534 

Ditto 

65 

80 

14.4 

0016 

960 

Ditto 

065 

24.9 

009 

1  08/1  05 

7.31/  7  52 

Ditto 

65 

*0 

14.4 

0017 

847 

Dnto 

1  30 

136 

0  40 

1.03,1  02 

12  85  12  98 

Diito 

65 

80 

144 

0017 

847 

Ditto 

2  50 

85 

1.40 

1  02/1  01 

17  29  17  46 

Two  valors  are  oomputed  for  and  tfS'acpj  corresponding  to  the  upper  and  lower  bounds  for  i  d  «  ar.  d  »  «v. 
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TABLE 


Data  summary  for  Cnchlow  und  Weils  for  titaniumHilloy  its  sets 


Material  properties 

Vessel  geometry4*' 

Test  conditions 

Calculated  values1*' 

Designation 

a,,  ksi 

ov.  ksi 

R,  in. 

t,  in. 

R/t 

Temp.  C 

c,  in. 

Oy.  ksi 

IR  so.) 

Vs 

lO  '-lin  ’HIh  *> 

Ti-8AC-1M0-!V 

138 

149 

15 

0.05 

300 

Room  Temp 

2.25 

30.0 

1.125 

1.06/1.05 

1.48/1.49 

Ditto 

138 

149 

15 

0.05 

300 

Ditto 

3.40 

25.0 

2-568 

1.08/106 

1.38/1.41 

Ditto 

138 

149 

1$ 

0.05 

300 

Ditto 

4.20 

19.0 

3.920 

1.06/1.05 

1  98/2.0 

Ditto 

.‘38 

149 

15 

0.05 

300 

Ditto 

8.0 

113 

14220 

1.07/1.06 

2.91/294 

Ditto 

138 

149 

15 

0.05 

300 

Ditto 

92 

10 

18.800 

1.07/1.06 

3.23/3.26 

Ditto 

138 

149 

333 

0.03 

1100 

Ditto 

16 

30 

0.30 

1.02/1.02 

1.33/1.33 

Ditto 

138 

149 

33.3 

0.03 

1100 

Ditto 

5.5 

15 

t.36 

1.02/1.02 

254/2.54 

Ditto 

138 

149 

70 

0.03 

2325 

Ditto 

4.25 

30 

0.18 

1.05/1X14 

0.79/8.0 

Ditto 

138 

149 

70 

0.03 

2325 

Ditto 

9.20 

15 

0.86 

1/1 

1.53 

'*'  The  R  =»  1 5-m  tests  involved  cylindrical  vessels,  the  R  -  33.3-in.  and  70-in.  tests  involved  curved  panels. 
Two  values  are  computed  for  and  of,1  nctpy  corresponding  to  the  upper  and  lower  bounds  for 

table 

DatJ  summary  for  Anderson  and  Sullivan 

for  titaniunt'alloy  vessels 

Material  properties 

Vessel  geometry 

Test  conditions 

Calculated  values'*’ 

Designation 

Oy,  ksi 

ov,  ksi 

R.  in. 

t,  in. 

*/■ 

Temp,  C 

c,  in. 

*2.  ksi 

50(c'/*')  llnh 
W»> 

Vs 

to’  «*,)-' 

10 

Ti-5Al-2.5  Sn 

193 

220 

30 

0.02 

150 

-196 

0.06 

190.4 

002 

sc/1.76 

0/0.83 

Ditto 

193 

220 

3.0 

002 

150 

-  196 

0.12 

164.9 

008 

1.90/1.49 

0.51/0.65 

Ditto 

193 

220 

3.0 

0.02 

150 

-1% 

0.11 

156.6 

006 

1.63/1.38 

072/0.85 

Ditto 

193 

220 

3.0 

0.02 

150 

-196 

022 

115.5 

026 

1.32/1.14 

0.82/0.95 

Ditto 

193 

220 

3.0 

0.02 

150 

-1% 

0.23 

105.1 

029 

1.26/1.18 

0.99/1.06 

Ditto 

193 

220 

3.0 

002 

150 

-196 

037 

84.9 

0.75 

1.25/1.18 

095/J.01 

Ditto 

193 

220 

3.0 

0.02 

150 

-196 

0.38 

74.6 

080 

U  8/1. 13 

1.27/1.33 

Ditto 

193 

220 

3.0 

002 

150 

-196 

0.47 

71.8 

1.22 

1.24/1.16 

1.05/ 1.1 3 

Ditto 

193 

220 

3.0 

0.02 

150 

-196 

0.49 

662 

1.33 

1.21/1.14 

1.22/1.30 

Ditto 

193 

220 

3.0 

0.02 

150 

-196 

0.74 

44.1 

3.00 

L 16/ 1.11 

1.90/1.99 

Ditto 

193 

220 

3.0 

0.02 

150 

-196 

0.73 

35.9 

2.90 

1.09/1.07 

3.10/3.16 

Ditto 

219 

240 

3.0 

0.02 

150 

-253 

0.M 

171.5 

001 

1.46/1.34 

1.61/1.75 

Ditto 

219 

240 

3.0 

0.02 

150 

-253 

007 

160.9 

003 

1.39/1.30 

1.14/1.22 

Ditto 

219 

240 

3.0 

002 

150 

-253 

0.09 

133.9 

005 

1.23/1.18 

1.51/1.38 

Ditto 

219 

240 

3.0 

0.02 

150 

-253 

0.13 

1214 

010 

1.20/1.16 

1.30/1.34 

Ditto 

219 

24 

3.0 

0.02 

150 

-253 

0.14 

1142 

0.10 

1.161. 13 

1. 50/1.54 

Ditto 

219 

240 

30 

002 

150 

-253 

0.26 

146 

037 

1  11/1.09 

1  54/1.56 

Ditto 

219 

240 

3.0 

0.02 

150 

-253 

0.24 

760 

032 

1.08/1.06 

2.12/2.15 

Ditto 

219 

240 

3.0 

0.02 

150 

-253 

040 

636 

088 

1.09/1.08 

1.79/1.82 

Ditto 

219 

240 

3.0 

0.02 

150 

-253 

0.47 

63.0 

122 

1  13/1  09 

1.52/1.55 

Ditto 

219 

240 

3.0 

002 

150 

-253 

0.38 

61.5 

080 

1.08/1.07 

105/206 

Ditto 

219 

240 

30 

002 

150 

-253 

0.49 

516 

1.33 

1.08/1.06 

2  25/230 

Ditto 

219 

240 

30 

0.02 

150 

-253 

080 

40.6 

355 

1.11/109 

217/221 

Ditto 

219 

240 

30 

0.02 

150 

-253 

0.78 

37.7 

338 

108/107 

265/268 

*•’  Two  values  are  computed  for  g>,  and  otfncvs  corresponding  to  the  upper  and  lower  hounds  for  » 
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TABLE 

Data  summary  for  Stchler  anJ  H'i Hums  for  brass  icsse/s 


Material  properties  Vessel  geometry  Test  conditions  Calculated  values'" 


Designatioa 

at,  ksi 

<*v>  ksi 

R,  in. 

f,  in. 

ft.l 

Temp,  C 

r,  in. 

»!• 

>G(<  •  ft1)  unh 
(ft  50f) 

•Pi 

i°H}  *C<P>)~  1 

10  10{m.,)(lb. 

Brass1** 

45 

57 

1.5 

0.00 1 

1500 

Room  Temp 

003 

50.6 

002 

x/1.86 

0/22.3 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

004 

454 

003 

*/l  52 

0/254 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

0.06 

37  1 

008 

1  74/1  71 

22  18/29  47 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

0.10 

30.0 

0  22 

1*44  1.21 

24.53;  29.2 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

010 

30.0 

022 

1.44,  1.21 

24  13/292 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

0.12 

255 

0  32 

1.31/1.16 

31  15/35.18 

Ditto 

45 

57 

1.5 

0  001 

1500 

Ditto 

0.12 

25.7 

0.32 

1  31/1.16 

30.65/34  62 

Ditto 

45* 

57 

1.5 

0.001 

1500 

Ditto 

013 

268 

0.37 

1  39,1  19 

24.55  24.68 

Ditto 

45 

57 

1.5 

0001 

1500 

Ditto 

0.15 

24.0 

0  50 

1.33,1.17 

27.74'31.53 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

015 

21.7 

0.50 

1.25,1.13 

36.03  39.86 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

0  17 

19.1 

0.64 

1.20  1.11 

42.73,46  2 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

0  21 

19.1 

0.98 

1.28/1  15 

32.55. 36  23 

Ditto 

45 

57 

1.5 

0001 

1500 

Ditto 

0  25 

14  2 

1.38 

1.17  1.10 

54  09  57.53 

Ditto 

45 

57 

1.5 

0.001 

1500 

Ditto 

0  30 

133 

200 

1.21 'Ml 

49  48  53.94 

Ditto 

45 

57 

1.5 

0001 

1500 

Ditto 

035 

11.0 

2.70 

1.16/1.09 

6481,6897 

Ditto 

40 

53  7 

15 

0.0U1 

2500 

Ditto 

0.032 

480 

0008 

x/1.86 

0  23  17 

Ditto 

40 

53.7 

25 

0.001 

2500 

Ditto 

0032 

46.1 

0.008 

r/1.66 

0  2S  14 

Ditto 

40 

5:  7 

25 

0.001 

2500 

Ditto 

0  062 

400 

003 

ac/1.40 

0  22  89 

Ditto 

40 

53.7 

2.5 

0.001 

2500 

Ditto 

0  094 

370 

0  07 

2.6  .1.36 

9.52  182 

Ditto 

40 

53.7 

15 

0.001 

2500 

Ditto 

0094 

38  7 

0.07 

ac/1.40 

0  16  16 

Ditto 

40 

53.7 

15 

0001 

2500 

Ditto 

010 

350 

008 

1.98;  1.30 

13  11  1998 

Ditto 

40 

53.7 

15 

0001 

2500 

Ditto 

0125 

32.2 

012 

1.76  1.26 

13  96  19  50 

Ditto 

40 

53.7 

15 

0.001 

2500 

Ditto 

0125 

30  2 

012 

1.54,1.21 

18.13  :3.C> 

Ditto 

40 

53.7 

15 

0.001 

2500 

Ditto 

0  15 

27.5 

0  18 

1.44.1.18 

19.50  2*  8  ) 

Ditto 

40 

53.7 

15 

0.001 

2500 

Ditto 

0156 

2fe2 

0  20 

1.39.1  16 

21  -I  25  65 

Ditto 

40 

53.7 

15 

0001 

2500 

Ditto 

0  156 

2*0 

0  2C 

I  43  l.|8 

19.5$  23  *3 

Ditto 

52 

632 

15 

0  003 

*33 

D.tto 

0048 

564 

oo;s 

x/l  86 

0  11  2 

Ditto 

52 

63.2 

25 

0.003 

833 

Dttto 

009$ 

500 

00*6 

r ,  1  60 

0  8  11 

Ditto 

52 

63  2 

15 

0.003 

833 

Ditto 

0  115 

48  4 

0.  !0 

x,l  54 

0  *  6  7 

Ditto 

52 

63.2 

15 

0.003 

833 

Ditto 

0  125 

42  5 

0  12 

x.,1  54 

0  9  1* 

Ditto 

52 

63.2 

15 

0003 

833 

Ditto 

0  125 

450 

0  12 

2  24  I  44 

5>f  8" 

Ditto 

52 

63.2 

15 

0003 

833 

Ditto 

0  125 

43  2 

0  »2 

l  RS  1  22 

8  25  11  18 

Ditto 

52 

63  2 

15 

0003 

833 

Ditto 

0  148 

41  3 

or 

1.85  1  38 

8  10  9  1? 

4*'  Two  talues  are  computed  for  <p}  and  oj,2  srv>»  corresponding  to  the  upper  ar.d  lower  bounds  for  c  .  6  *  cr.  c  rt 
,M  Brass  shim  stock. 
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